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SUMMARY 


This  report  presents  results  of  measurements  and  calculations  on  backscat- 
tering  of  intermediate-  and  high-frequency  ultrasonic  waves  by  42  voids 
and  inclusions  embedded  in  titanium  and  1N100  alloys.  The  primary  program 
objective  was  to  evaluate  the  Born  inversion  technique  (BIT)  in  the  forms 
developed  at  the  Ames  Laboratory,  Iowa  State  University,  and  the  Rockwell 
Science  Center  against  nternal  bulk  flaws  under  simulated  field  conditions 

The  program  consisted  of  five  tasks.  The  first  task  was  the  adaptation  of 
the  Born  inversion  software  and  testing  of  the  assembled  hardware  system 
using  the  adapted  software  on  flaws  with  known  characteristics.  The  second 
task  involved  the  fabrication  by  Pratt  &  Whitney  Aircraft  Group  of  titanium 
alloy  and  powder  metallurgy  nickel-based  superalloy  test  specimens  contain¬ 
ing  40  volumetric  flaws  in  the  250  to  1,300  um  size  range.  The  third  task 
included  the  development  and  qualification  of  two  sets  of  flaw  character¬ 
ization  procedures  on  the  test  flaws.  The  BIT  was  incorporated  into  the 
primary  flaw  characterization  procedures  along  with  the  spectrum  splicing 
technique  (SST)  that  combines  the  backscattered  amplitude  data  from  the 
medium-  and  high-frequency  regions  to  obtain  quickly  an  indication  of  the 
appropriateness  of  a  transducer  for  acquiring  the  waveforms  for  Born  inver¬ 
sion.  Prior  to  any  metallographic  sectioning,  the  test  flaws  were  also 
characterized  as  to  size,  shape,  and  composition  using  procedures  that 
incorporated  the  phase-comparison  technique  and  the  satellite-pulse  obser¬ 
vation  technique  (SPOT).  The  final  two  tasks  involved  metallographic  sec¬ 
tioning  of  three  real  manufacturing  flaws  and  comparison  of  the  BIT  and 
SPOT  size  estimates  with  the  nominal /actual  flaw  dimensions. 

It  is  concluded  that  multiple  measurements  made  with  waves  of  different 
mode,  incidence  angle,  and  dominant  frequency  and  no  less  than  three  phys¬ 
ical  models  for  flaw-ultrasound  interactions  (one  for  voids  and  the  other 
two  for  weakly  and  strongly  scattering  inclusions)  are  necessary  to  solve 
the  inverse  problem  for  flaw  type  and  size  in  an  unambiguous  manner. 
Experiments  with  known  flaws  guided  the  model  development  efforts.  The 
Born  diameter  estimates  were  determined  to  be  independent  of  flaw  composi¬ 
tion,  size,  and  shape,  but  linearly  related  to  the  period  of  the  ultrasonic 
wave.  The  BIT  yielded  size  estimates  with  barely  acceptable  accuracy  for 
those  voids  and  inclusions  for  which  a  transducer  with  a  center  frequency 
near  the  periodicity  of  the  flaw's  magnitude  spectrum  was  available  (coin¬ 
cidence  rule).  The  BIT  completely  broke  down,  however,  for  22  of  the  flaws 
In  the  IN100  test  specimens  because  these  flaws  could  not  be  detected  by 
more  than  one  transducer.  The  coincidence  rule  was  violated  for  most  of 
these  flaws. 

External  satellite  pulses  originating  from  the  waves  that  circumnavigated 
the  flaws  were  observed  for  each  flaw.  Additional  (internal)  satellite 
pulses,  on  the  other  hand,  were  observed  only  for  the  tungsten  carbide 
inclusions  in  the  titanium  alloy  specimens.  The  type  and  size  of  the 
flaws  were  determined  by  comparing  the  flaw  "signatures"  (time-domain 
responses)  and  flaw  "profiles"  (frequency-domain  responses)  to  model 


predictions*  The  sizing  errors  for  the  36  detectable  voids  and  inclusions 
in  the  diffusion-bonded  specimens  were  such  that  merely  100  um  (4  mils) 
need  to  be  added  to  the  SPOT  estimate  in  order  to  be  almost  certain  (i.e*, 
95  percent  sure)  that  the  actual  flaw  size  does  not  exceed  this  conserva¬ 
tive  estimate*  The  SPOT  size  estimates  were  bracketed  by  the  maximum 
linear  dimensions  of  the  three  metallographically  sectioned  real  flaws. 
Semiautomated  satellite  analysis  procedures  (SSAP)  that  waive  the  require¬ 
ment  for  a  highly  trained  ultrasonic  examiner  were  developed  and  qualified 
on  the  32  inclusions. 
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1 .  INTRODUCTION 


1 . I  Background 


Turbine  engine  component  life-extension  technology  now  being  devel¬ 
oped  under  the  name  "retirement  for  cause"  requires  reliable  and  accurate 
ultrasonic  techniques  for  characterizing  material  flaws  as  to  type,  size, 
location,  shape,  orientation,  and  composition.  Ultrasonic  techniques  based 
solely  on  signal  amplitude  (i.e.,  the  amplitude-comparison  and  amplitude- 
drop  techniques)  are  incapable  of  characterizing  flaws  with  the  accuracy 
needed  for  component  life  prediction  systems  based  on  fracture  mechanics. 
The  limited  capability  of  present-day  ultrasonic  techniques  has  become  a 
barrier  to  the  fulfillment  of  inspection  requirements  in  a  number  of  situa¬ 
tions  of  interest  to  the  Air  Force.  In  situations  where  detailed  informa¬ 
tion  about  a  detected  flaw  could  be  used  in  component  life  extension,  a 
Quantitative  Ultrasonic  Flaw  Characterization  Module  would  be  a  necessary 
supplement  to  the  qualitative  capability  now  available. 

One  objective  of  the  Interdisciplinary  Program  for  Quantitative  Non¬ 
destructive  Evaluation,  established  jointly  by  the  Air  Force  Materials  Lab¬ 
oratory  (AFML)  and  the  Defense  Advanced  Research  Projects  Agency  (DARPA)  in 
the  mid-1970's,  was  to  develop  quantitative  inspection  techniques.  Ele¬ 
ments  of  this  extensive  development  effort  became  available  by  1980  for 
transition  into  production  and  depot  ultrasonic  inspection  equipment.  Two 
tasks  remained,  however,  prior  to  reducing  to  practice  the  most  mature  of 
the  inverse  scattering  algorithms  developed  under  earlier  AFML/ DARPA  con¬ 
tracts.  The  first  task  was  to  develop  procedures  featuring  the  Born  inver¬ 
sion  technique  (BIT)  supplemented  by  a  transducer  selection  protocol.  The 
second  task  was  to  qualify  the  Born  inversion  procedures  (BIP)  in  the  forms 
developed  at  Ames  Laboratory  (1^)  and  Rockwell  Science  Center  (2)  against  a 
statistically  large  set  of  simulated  and  real  volumetric  flaws  (voids  and 
inclusions).  The  Test  Bed  Program  (2^)  resulted  in  a  set  of  Born  inversion 
procedures  that  requires  a  preliminary  look  at  the  modulated  magnitude 
spectrum  of  the  flaw  H(f)  to  insure  that  the  bandwidth  of  the  transducer 
is  centered  on  the  portion  of  H(f)  containing  the  first  peak.  These  proce¬ 
dures  were  successfully  applied  to  a  natural  flaw  in  a  powder  metallurgy 
nickel-based  superalloy  sample.  In  1981,  Southwest  Research  Institute 
(SwRI)  initiated  the  Flaw  Characterization  Module  Development  Program  to 
further  test  the  Born  inversion  procedures  against  a  large  set  of  flaws 
under  simulated  field  conditions. 

1.2  Objectives 


The  first  objective  of  the  Flaw  Characterization  Module  Program  was 
to  test  the  supplemented  BIT  on  a  large  number  of  representative  flaw  types 
(not  less  than  40).  Additionally,  the  approximate  dimensions  of  the  flaws 
were  to  be  determined  by  at  least  one  other  nondestructive  means  ("second 
method")  prior  to  any  metallographic  sectioning.  The  second  objective  was 
to  implement  the  BIP,  if  proven  valid,  in  a  module  compatible  with  state- 
of-the-art  ultrasonic  equipment. 


The  first  objective  of  the  program  was  accomplished  in  five  tasks. 

A  flow  diagram  of  the  work  performed  to  qualify  the  BIP  for  simulated  and 
real  internal  bulk  flaws  is  shown  in  Figure  1.  Task  1  consisted  of  data 
acquisition,  storage,  and  processing  system  ("hardware")  assembly;  zero- 
of-time  algorithm,  Born  inversion  algorithm,  and  transducer  selection  pro¬ 
tocol  ("software")  adaptation;  and  functional  testing  of  the  hardware- 
software  combination  using  spherical  voids  of  known  size  in  standard  speci¬ 
mens.  Task  2  involved  the  fabrication  of  test  specimens  of  typical  air¬ 
craft  engine  materials  containing  at  least  40  isolated  bulk  flaws  with 
characteristics  unknown  to  the  ultrasonic  examiners.  Task  3  consisted  of 
adaptation  and  further  development  of  the  BIP;  testing  of  these  procedures 
against  the  flaws  fabricated  or  collected  during  Task  2;  and  application 
of  additional  sizing  techniques.  After  all  ultrasonic  examinations  of  the 
flaws  in  the  test  specimens  were  completed,  half  of  the  collected  manufac¬ 
turing  (real)  flaws  were  metallographically  sectioned  to  determine  the 
actual  flaw  characteristics  (Task  4).  Task  5  involved  the  comparison  of 
the  Born  diameters  and  other  size  estimates  to  the  nominal  diameters  of  the 
flaws  in  the  diffusion-bonded  specimens  and  the  actual  dimensions  of  the 
sectioned  real  flaws  detected  in  rejected  parts.  The  second  program  objec¬ 
tive  was  not  accomplished  because  the  Born  diameter  results  indicated  that 
it  would  be  premature  to  construct  a  quantitative  flaw  characterization 
module . 


1.4  Report  Organization 

The  remaining  part  of  the  report  summarizing  the  work  accomplished 
during  the  program  tasks  and  subtasks  is  divided  into  seven  sections. 
Table  1  shows  the  organization  of  the  report  along  the  specific  tasks  and 
subtasks  for  ready  reference. 

Table  1 

ORGANIZATION  OF  REPORT  ALONG  PROGRAM  TASKS  AND  SUBTASKS 


Task/Subtask  and  Title 


1.1 

Technique 

>  Review 

1.2 

Software 

Transfer 

1.3 

Hardware 

Assembly 

1.4 

Software 

Adaptation 

1.5 

Protocol 

Development 

2. 

Specimen 

Fabrication 

3. 

Protocol 

Testing 

4. 

Metallographic  Sectioning 

5. 

Analysis 

of  Results 

Report  Section  and  Title 

2.  Characterization  Techniques 

3.  Software  Transfer  and  Adapation 

4.  Procedure  Development 

5.  Test  Specimens 

6.  Procedure  Qualification 

7.  Discussion 


8.  Summary 
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STANDARD  SPECIMENS 
Task  1 


BACKSCATTERED 

WAVEFORMS 


Tasks  1  &  3 


TEST  SPECIMENS 


Task  2 


End  of 


Task  3 


AMPLITUDE  COMP.  TECH. 
PHASE  COMP.  TECH. 
AMPLITUDE  DROP  TECH. 
SPECTRUM  SPLICING  TECH. 


Task  4 


METALLOGRAPHIC 

SECTIONING 


ACTUAL  FLAW 
CHARACTERISTICS 


FLAW  COMPOSITION 
PROPER  TRANSDUCERS 


BORN  INVERSION 
TECHNIQUE 


BORN 

1 

DIAMETERS 

h 

ADDITIONAL  SIZING 
TECHNIQUES 


SIZE 

ESTIMATES 


COMPARISON 


Task  5 


CONCLUSIONS 


Figure  1.  Flow  Diagram  of  Effort  to  Qualify  the  Supplemented  Born 
Inversion  Technique  (BIT)  for  Voids  and  Inclusions 


2.  FLAW  CHARACTERIZATION  TECHNIQUES 


2.1  General 


The  inverse-scattering  problem  is  an  ill-posed  mathematical  prob¬ 
lem  (3 ,h) .  That  is,  flaws  that  differ  significantly  in  size,  shape,  com¬ 
position,  etc.,  may  produce  backscattered  waveforms  y(t)  that  are  essen¬ 
tially  the  same.  To  obtain  reliable  identification  of  flaw  type  and  accu¬ 
rate  estimates  of  flaw  size,  shape,  orientation,  and  composition,  this 
weakness  of  the  inverse  problem  must  be  recognized  and  appropriately  dealt 
with.  Quantitative  flaw  definition  requires  a  good  understanding  of  the 
interaction  of  ultrasound  with  flaws  of  various  characteristics.  This 
interaction  is  modeled  in  Figure  2  for  an  inclusion-like  (penetrable)  flaw 
with  smooth  boundaries.  The  biggest  challenge  for  inverse-scattering 
approaches  to  flaw  characterization  is  the  establishment  of  unambiguous 
cause-and-effect  relationships  among  the  fracture-related  flaw  characteris¬ 
tics  (X^,  i=l,  2,...M)  and  the  various  extractable  signal  parameters  (Y^, 
J-l,  2,...N). 

It  is  well  known  that  both  the  frequency  and  angular  distributions 
of  the  backscattered  wave  amplitudes  yn(?,t)  depend  on  the  characteristics 
of  the  flaw.  To  solve  the  inverse  problem,  the  distribution  of  the  ampli¬ 
tudes  of  the  backscattered  waves  must  then  be  measured  over  a  set  of  wave 
modes  (longitudinal  vs.  shear),  transducer  center  frequencies  (medium  vs. 
high),  and  beam  angles.  It  follows  that  the  inversion  process  sought  is 
one  that  allows  estimation  of  the  ith  flaw  characteristic  on  the  basis  of 
just  a  few  independent  measurements  (i.e.,  signal  parameters). 

If  one  knows  a  priori  that  the  flaw  is  a  spherical  void,  only  one 
signal  parameter  needs  to  be  extracted  from  y(t)  to  estimate  the  one 
unknown  flaw  characteristic  (diameter).  The  number  of  flaw  characteris¬ 
tics  for  a  nonspherical  inclusion,  on  the  other  hand,  is  no  less  than  five. 
When  M  is  large,  determination  of  the  effects  of  each  on  a  given  signal 
parameter  Yj  requires  a  combination  of  analytical,  experimental,  and  stat¬ 
istical  approaches.  One  cause  (X^)  may  have  several  effects  (Yj),  and  the 
same  effect  may  be  due  to  several  causes. 

The  BIT  is  a  method  for  estimating  the  geometrical  characteristics 
(shape,  size,  and  orientation)  of  internal  bulk  flaws  and  uses  as  input 
data  the  specularly  reflected  wave  [ y ^ ( * )  in  Figure  2(a)].  This  technique 
needs  supplementation  not  only  in  establishing  the  flaw's  type  and  composi¬ 
tion,  but  also  in  making  sure  that  the  transducer's  bandwidth  requirement 
(Born  bandwidth,  0.5<ka<2.5)*  is  met.  The  three  amplitude  and  phase  analy¬ 
sis  techniques  supplementing  the  BIT  are  described  first.  This  is  followed 
by  discussions  of  the  BIP  (including  a  transducer  selection  protocol)  and 
an  additional  sizing  technique  that,  unlike  the  BIT ,  utilizes  the  trailing 
signals  ("satellite  pulses")  resulting  from  the  secondary  backscattered 
waves  [ y2 ( * )  and  y^(*)  in  Figure  2(a)]  in  the  estimation  of  flaw  size. 

*k(“2 n/x)  is  wave  number,  X(»c/f)  is  wavelength,  c  is  sound  velocity  of 
host  material,  f  is  transducer  center  frequency,  and  a  is  flaw  radius. 


y, (r\t)  Reflected 

yiir.t)  Internal 

y,(r*,t)  External  (Creeping) 


Backscattered  Waves 

(a)  Physical  Model  of  Flaw-Ultrasound  Interaction 


FLAW  WAVEFORMS 

y(t)  _  w(t)  Reference 

^  y(t)  Flaw 

* 


▼ 

Causes  -  Effects  - 

Flaw  Characteristics  Signal  Parameters  Y- 

i  -  1,  2....M  j  =  l,  2....N 

(b)  Linear  System  Model  of  Flaw-Ultrasound  Interaction 

Figure  2.  Models  for  Predicting  the  Response  of  a  Penetrable 
Volumetric  Flaw  to  Ultrasound 


2.2  Amplitude-Comparison  Technique 

Analysis  of  the  maximum  amplitude  of  the  rectified  flaw  waveform 
jy(t) |  is  most  useful  for  flaw  detection.  The  size  of  flaws  that  are  small 
compared  to  the  wavelength  (ka<0.5,  Rayleigh  region)  can  be  estimated  by 
the  amplitude-comparison  technique.  However,  even  in  the  Rayleigh  region, 
the  maximum  signal  amplitude  depends  strongly  on  flaw  composition.  Only 
in  very  favorable  circumstances  will  the  maximum  signal  amplitude  bear  an 
unambiguous  relationship  to  flaw  size. 

The  amplitude-comparison  technique  is  more  reliable  for  distinguish¬ 
ing  among  voids  and  inclusions  than  for  flaw  sizing.  It  is  known  that, 
due  to  their  internal  waves,  tungsten  carbide  and  other  weakly  scattering 
inclusions  exhibit  consistently  more  structured  "signatures"  (time-domain 
responses)  and  more  modulated  "profiles"  (frequency-domain  responses)  than 
voids  (5,6^).  For  a  given  void-transducer  combination,  the  creeping-wave- 
to-ref lected-wave  amplitude  ratio  improves  in  going  from  longitudinal  to 
shear  waves  (7) . 

2.3  Phase-Comparison  Technique 

This  technique  provides  a  more  quantitative  means  than  the  amplitude- 
comparison  technique  of  determining  the  flaw  type  and  composition.  Inclu¬ 
sions  with  a  lower  acoustic  impedance  Z*  than  that  of  the  host  material — 
void  is  a  special  case  of  this  type  of  inclusion — can  be  readily  distin¬ 
guished  from  those  inclusions  that  have  a  higher  acoustic  impedance  than 
the  host  material  by  observing  the  phase  of  the  specularly  reflected  pulse 
["locomotive"  of  the  pulse  train  (4,8,9,10)] .  Additional  flaw-composition 
information  can  be  extracted  from  the  signs  of  the  satellite  pulses 
("trailers")  accompanying  the  reflected  pulse.  This  will  be  further  dis¬ 
cussed  in  Subsection  2.7. 

2 . 4  Amplitude-Drop  Technique 

This  technique  is  commonly  used  in  the  field  for  flaw  sizing.  When 
the  largest  dimension  of  an  elongated,  needle-like  flaw  exceeds  that  of 
the  ultrasonic  beam  incident  perpendicularly  on  it,  there  are  generally 
not  many  sizing  problems.  The  length  of  a  4-mm  long  natural  Inclusion  was 
accurately  estimated  by  the  amplitude-drop  technique  using  a  focused  trans¬ 
ducer  ( 2 _) .  The  through-wall  dimension  of  this  flaw  (about  140  ;.m)  was, 
however,  overestimated  by  a  factor  of  5  by  the  focused  C-scan  technique. 

The  correlations  were  also  poor  for  the  size  estimates  obtained  by  the 
focused  C-scan  technique  for  the  small  inclusions  (d<^500  urn)  in  some 
titanium  powder  compacts  (9). 

2 . 5  Born  Inversion  Procedures 

The  BIT,  credited  to  Rose  and  Krumhansl  (H_),  requires  a  single 
waveform  from  an  appropriate  transducer  to  obtain  a  size  estimate 


oc  where  p  is  the  density  of  the  inclusion 
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(Born  diameter  D)  corresponding  to  the  projection  of  the  flaw  along  the 
direction  of  insonif ication.  The  transducer  selection  protocol  features 
the  spectrum  splicing  technique-  The  shape  and  orientation  of  the  flaw 
may  be  estimated  by  moving  the  pulse-echo  transducer  around  the  flaw  and 
inverting  the  waveforms  received  from  several  directions  or  by  using  a 
multi-viewing  transducer  (12) . 

The  input  to  the  Born  inversion  algorithm  (BIA)  consists  of  the  pro¬ 
perly  time-shifted  real  part  of  the  deconvolved  flaw  spectrum  [A(f)].  A 
deconvolution  algorithm  of  the  Wiener-filter  type  was  derived  to  minimize 
the  effects  of  transducer  and  material  properties  on  the  output  of  the  BIA 
(i.e.,  the  Born  diameter).  A  zero-of-time  algorithm  (ZOTA)  was  adapted 
for  locating  the  center  of  the  flaw  along  the  observation  direction.  The 
spectrum  splicing  technique  and  the  three  algorithms  incorporated  into  the 
BIP  are  described  in  the  following  paragraphs. 

2.5.1  Transducer  Selection  Protocol 

The  most  serious  limitation  of  the  BIT  is  the  requirement 
for  a  wideband  ultrasonic  transducer  that  is  effective  over  the  specific 
0.5  to  2.5  ka  range.  This  bandwidth  requirement  is  at  the  limit  of  state- 
of-the-art  transducer  technology.  The  limitation  on  the  bandwidth  dic¬ 
tates  the  center  frequency  of  the  "proper"  transducer  (Born  frequency)  to 
be 

„  c(0.5  x  2. 5)**  _  2,200  (I) 

TTd  "  d 


where  d  is  flaw  diameter  in  um,  c  =  6.1  km/sec  for  the  titanium  and  IN100 
alloys  of  interest,  and  F  is  in  MHz.  The  transducer  selection  problem  is 
solved  by  examining  the  deconvolved  flaw  spectrum  [magnitude  spectrum  H(f)] 
obtained  by  using  two  or  three  longitudinal-wave  transducers  and  splicing 
the  individual  spectra  together  (2). 

The  spliced  spectrum  expected  for  a  1,000-um  diameter  void  is 
shown  in  Figure  3  and  that  expected  for  a  500-um  diameter  void  is  shown  in 
Figure  4.  Note  that  the  modulations  of  the  theoretical  magnitude  spectra 
are  quite  regular  (10)  and  the  frequency  of  the  first  peak  is  nearly  equal 
to  p,  the  spectrum  periodicity.  Note  also  that  to  a  good  approximation 
the  product  of  p  and  d  is  a  constant  (2,200).  It  follows  that  the  Born 
frequency  for  a  strongly  scattering  ("strong")  flaw  is  given  by  the  aver¬ 
age  periodicity  of  the  composite  magnitude  spectrum  obtained  by  splicing 
together  the  individual  spectra  of  two  or  three  intermediate-  and  high- 
frequency  transducers.  This  is  the  essence  of  the  spectrum  splicing  tech¬ 
nique  (SST).  The  bandwidth  of  the  Born  transducer  not  only  must  be  suit¬ 
ably  large,  but  must  also  be  centered  on  the  first  peak  of  H(f).  It 
follows  from  the  generalization  of  this  coincidence  rule  that  the  Born 
frequency  for  a  weakly  scattering  ("weak”)  flaw  is  given  by  70  percent  of 
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the  periodicity  of  its  spliced  magnitude  spectrum.*  In  summary,  the  Born 
frequency  is  given  by 


F  = 


p  for  strong  scatterers 


(  0.7p  for  weak  scatterers 

where  p  is  the  periodicity  of  the  flaw's  spliced  magnitude  spectrum. 
2.5.2  Deconvolution  Algorithm 


(2) 


Deconvolution  to  remove  the  effects  of  the  data  acquisition 
system  from  the  flaw  waveform  y(t)  is  the  first  preprocessing  step  in  the 
application  of  the  BIT.  All  deconvolution  algorithms  are  sensitive  to 
noise.  The  deconvolution  process  can  be  performed  in  the  frequency  or  the 
time  domain.  In  a  recently  developed  time-domain  deconvolution  algo¬ 
rithm  (4),  a  trade-off  is  made  between  the  ability  to  detect  the  trailing 
signals  associated  with  the  specularly  reflected  pulse  (resolution)  and 
insensitivity  to  material,  electronic,  and  other  types  of  noise.  A  Wiener- 
filter  type  deconvolution  that  is  performed  in  the  frequency  domain  is 
described.  It  follows  from  the  system  model  of  Figure  2(b)  that,  in  the 
absence  of  noise,  the  flaw  spectrum  Y(f)  is  simply  equal  to  the  product  of 
the  flaw's  transfer  function  H(f)  and  the  reference  spectrum  W(f).  The 
front  or  back  surface  of  the  test  specimen  containing  the  flaw  is  an 
acceptable  reference  reflector.  A  Wiener-filter  type  deconvolution  algo¬ 
rithm  of  the  form 


H  (f )  = 


Y(f)  W*(f) 
|W(f)|Z  +  Q2 


(3) 


where  the  asterisk  denotes  complex  conjugate  and  Q2  is  the  average  noise 
power  density,  can  desensitize  the  results  to  the  presence  of  noise  at 
high  and  low  frequencies  where  the  transducer  is  inefficient.  Pulse  reso¬ 
lution  can  be  traded  for  Insensitivity  to  noise  and  vice  versa  in  the 
selection  of  the  Q  parameter. 


2.5*3  Zero-of-Tlme  Algorithm 

The  deconvolved  flaw  spectrum  that  results  from  performing 
the  operations  required  by  Equation  (3)  is  generally  complex.  The  phase 
error  2iTfx  must  be  eliminated  from  H(f)  before  it  can  be  input  to  the  BIA. 
This  requires  shifting  the  impulse  response  h(t)  by  x  so  that  the  centroid 
of  the  flaw  is  at  t  =  0.  If  the  origin  of  time  coincided  with  the  time  at 
which  the  fictitious  pulse  labeled  C  in  Figure  5  returned  from  the  center 
of  the  flaw  (point  Q0),  there  would  be  no  need  for  time  shifting  of  h(t). 
If  sufficient  low-frequency  (ka <  0.5)  information  was  available,  the 
required  time  shift  could  be  determined  by  extrapolating  the  phase 


*The  70  percent  reflects  the  differences  in  velocities  and  travel  paths 
between  the  externally  scattered  waves  of  a  strong  flaw  and  the  internally 
scattered  waves  of  a  weak  flaw  [see  Figure  2(a)  for  the  travel  path 
d  if  f erences] . 
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FLAW  CENTER 


Figure  5.  Physical  Model  for  the  Born  Inversion  Technique. 

The  reflected  pulse  (R)  is  real  and  the  center 
pulse  (C)  is  fictitious. 


spectrum  to  zero  frequency.  However,  state-of-the-art  Born  transducers 
(i.e.,  transducers  centered  at  ka  =  1.1)  do  not  permit  this  extrapolation 
process . 

The  method  recommended  by  Ames  Laboratory  (J^)  was  used  in  this  pro¬ 
gram  to  determine  the  zero  of  time  for  h(t).  This  method  is  based  on  the 
premise  that  if  H(f)  contains  a  phase  error  of  2irfT  then  the  peak,  of  the 
cross-sectional  area  function  for  the  flaw  a(z)  is  shifted  from  the  proper 
origin  by  cx/2.  Figure  6  shows  the  cross-sectional  area  function  and  the 
time-shifted  Impulse  response  obtained  with  a  "proper"  (F  =  6  MHz)  trans¬ 
ducer  for  a  400  x  800  urn  oblate  spheroidal  void  in  a  titanium  6-4  specimen 
(13).  Note  that  the  peak  of  a(z)  and  the  peak  of  h(t)  corresponding  to 
the  C  pulse  (second  half  cycle)  occur  at  the  same  location  (i.e.,  at  the 
centroid  of  the  impulse  response)  and  have  the  same  polarity.  The  first 
half  cycle  of  h(t)  has  the  opposite  polarity  and  is  identified  to  be  the  R 
pulse  of  Figure  5.  The  BIT  reduces  to  a  mere  time-of-f light  reconstruction 
for  the  R  pulse  since  it  does  not  utilize  the  arrival  time  for  the  size- 
information-carrying  creeping  wave  (external  satellite  pulse  Se)  clearly 
seen  in  both  y(t)  and  h(t)  of  Figure  6- 

2.5.4  Born  Inversion  Algorithm 

The  one-dimensional  BIA  produces  a  characteristic  function 
J(r)  from  the  real  part  of  the  flaw's  transfer  function  A(k)  according  to 
the  inverse  transformation 


J(r) 


constan 


A(k) 


sin  2kr 
2kr 


dk 


(4) 


The  characteristic  function  shown  in  Figure  6  has  zero  slope  at  the  origin. 
An  estimate  of  the  flaw  radius  (half  of  Born  diameter  D)  can  be  obtained 
from  the  characteristic  function  as  the  range  at  which  J(r)  drops  to 


one  half  of  Its  peak  value*  The  Born  diameter  can  also  be  calculated 
according  to  the  formula 


D  =  2  Area  under  J(r)  (5) 

Peak  value  of  J(r) 

Note  that  the  product  of  the  transducer  center  frequency  (6  MHz)  and  the 
Born  diameter  (360  urn)  is  about  2,200  mHz.  This  agrees  with  the  constant 
of  Equation  (1). 


Satellite-Pulse  Observation  Technique 


Volumetric  material  flaws  may  be  broadly  classified  as  penetrable  or 
impenetrable  (i.e.,  voids).  Penetrable  flaws  may  be  divided  into  weak 
flaws  (Born-like  inclusions)  and  strong  flaws  [Franz-like  inclusions  (14)} 
The  sketches  on  the  left  side  of  Figure  7  reveal  the  main  differences  in 
the  ultrasonic  scattering  mechanisms  of  voids,  weak  inclusions,  and  strong 
inclusions.  All  types  of  volumetric  flaws  can,  in  principle,  support  tan¬ 
gentially  radiating  Franz-type  "creeping"  waves  as  indicated  by  the  wiggly 
lines  between  the  shadow  boundary  points  and  Q2  of  the  various  flaws. 

In  addition  to  the  externally  scattered  (Se)  waves,  weak  flaws  can  also 
support  internally  scattered  (S^)  waves  resulting  in  a  pair  of  satellite 
pulses.  The  satellite-pulse  observation  technique  (SPOT)  is  based  on  the 
interpretation,  in  terms  of  flaw  types  and  dimensions,  of  the  separations 
a  in  time  of  arrival  between  the  specularly  reflected  R  pulse  and  its 
generally  ignored  trailers  (satellites  Sj  and  Se)  contained  in  the  flaw 
waveform  (15) .  The  analytical  and  supporting  experimental  SPOT  results 
obtained  for  the  three  flaw  types  of  interest  are  presented. 


2.6.1  Physical  Model  for  Voids 


It  is  known  that  the  small  pulses  launched  in  opposite  direc¬ 
tions  (i.e.,  clockwise  and  counterclockwise)  from  the  shadow  boundaries  of 
a  strong  scatterer  of  convex  shape  travel  in  unison  on  the  instrument 
screen  with  the  specularly  reflected  pulse  (hence  'satellite"  pulses). 

The  upper  diagrams  in  Figure  7  show  the  interaction  of  a  high-frequency 
longitudinal  or  shear  wave  with  a  void  resulting  in  a  reflected  pulse  (R) 
and  an  externally  scattered  pulse  (Se).  For  clarity,  only  the  clockwise 
creeping  wave  is  shown.  Each  time  this  creeping  wave  reaches  point  Q2, 
tangential  reradiation  in  the  form  of  a  tangentially  scattered  wave  from 
the  void  during  its  circumvention  is  launched  in  the  direction  to  be 
received  by  the  pulse-echo  transducer.  The  separation  of  the  first  satel¬ 
lite  from  the  R  pulse  is  given  by 


where  v  is  the  creeping-wave  velocity.  For  longitudinal -wave  (LW)  inci¬ 
dence,  c  »  6.1  km/sec  and  v  is  taken  to  be  0.87c  for  an  average-sized 
(800-i.m  diameter  standard)  void  insonified  with  a  standard  high-frequency 
(5-MHz)  transducer  (^_6 ) .  Published  creeping-wave  velocities  range  from 
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Figure  7.  Physical  Models  for  the  Interaction  of  High-Frequency 
Ultrasound  With  Volumetric  Flaws  Resulting  in  External 
and/or  Internal  Scattered  Waves  and  Satellite  Pulses 
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0.80c  to  0.95c  (_10,  n>  The  diameter  of  void-like  flaws  can  be  estimated 
on  the  basis  of  longitudinal-wave  o  measurements  by 

d  =  2.18  ae  (7) 

where  d  is  in  micrometers  and  oe  is  in  nanoseconds. 

For  shear-wave  (SW)  incidence,  c  =  3.1  km/sec  and  v  is  equal 
to  the  Rayleigh  wave  velocity  calculated  to  be  0.93c  (15) .  On  the  basis 
of  shear-wave  o  measurements,  the  void  diameter  is  estimated  by 

d  =  1.15  °e  (8) 

where  d  is  in  micrometers  and  oe  is  in  nanoseconds. 

The  composite  waveforms  obtained  by  interrogating  a  void  of 
unknown  size  and  shape  (i.e.,  one  of  the  test  flaws)  with  high-frequency 
longitudinal  and  shear  waves  are  reproduced  in  the  upper  diagrams  of  Fig¬ 
ure  8.  As  expected,  the  separation  of  the  SW  satellite  from  the  R  pulse 
is  nearly  twice  as  large  as  that  of  the  associated  LW  pulses.  Note  also 
that  the  reflected  and  the  externally  scattered  pulses  are  of  the  same 
polarity.  Both  the  leading  and  trailing  pulses  are  negative  because  there 
is  a  drop  in  the  acoustic  impedance  at  the  metal-void  boundary. 

2.6.2  Physical  Model  for  Tungsten  Carbide  Inclusions 

The  waveform  is  much  more  structured  when  both  externally 
and  internally  scattered  waves  are  produced  by  an  inclusion  with  sound 
velocities  nearly  the  same  as  those  of  the  host  alloy  matrix.  It  is  well 
known  that  the  back  side  of  a  tungsten  carbide  (WC)  or  other  Born-like 
inclusion  in  a  titanium  alloy  can  act  as  an  acoustic  lens  (10,18, 19) ♦ 
Because  of  this  focusing  effect,  the  amplitude  of  the  internally  scattered 
pulse  (SjL  in  the  middle  diagrams  of  Figures  7  and  8)  can  be  comparable  to 
that  of  the  R  pulse.  The  incident  longitudinal  wave  (I  in  Figure  7)  can 
propagate  through  the  middle  of  a  Born-like  inclusion  and  be  scattered 
from  its  back  surface  either  as  a  longitudinal  or  a  shear  wave  (18,20). 

The  same  mode-conversion  phenomena  can  also  be  observed  with  incident 
shear  waves  (_5  ) .  For  longitudinal-wave  incidence,  the  separation  of  the 
first  satellite  pulse  (S^)  from  the  R  pulse  is  given  by  2d/c  where 
c  =  6.6  km/sec  for  a  tungsten  carbide  inclusion.  The  diameter  of  a  WC 
inclusion  can  be  estimated  on  the  basis  of  longitudinal-wave  o  measure¬ 
ments  by 


d  =  3.30  o1  (9) 

where  d  is  in  micrometers  and  oj  is  in  nanoseconds. 

It  is  known  that  longitudinal  creeping  waves  can  propagate 
not  only  on  the  surfaces  of  strong  inclusions  but  also  on  the  surfaces  of 
weak  inclusions  (6^,21, 2_2).  If  observed,  the  separation  of  the  second 
satellite  pulse  (Se)  from  the  R  pulse  may  also  be  used,  in  accordance  with 
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Figure  8.  Analytical  Models  for  the  Estimation  of  Flaw  Size  d  in 
Terms  of  Doublet  Separations  oe  and  ot  Extracted  From 
the  Signatures  of  Three  Test  Flaws 


Equation  7,  to  estimate  the  size  of  a  WC  inclusion.  Because  the  density 
of  the  WC  inclusion  is  larger  than  that  of  the  titanium  alloy  matrix,  the 
R  pulse  is  expected  to  have  a  positive  polarity.  The  pulse  has  the 
same  polarity  and  the  Se  pulse  has  the  usual  negative  polarity. 

Externally  scattered  (creeping)  waves  appear  to  be  less 
important  when  the  WC  inclusion  is  interrogated  with  a  shear  wave .  The 
handbook  value  for  the  shear-wave  velocity  of  WC  is  4  km/sec.  However,  in 
one  study  of  the  scattering  behavior  of  WC  inclusions  (10)  the  shear  wave 
velocity  was  allowed  to  vary  between  2  and  4  km/sec.  We  decided  to 
determine  the  effective  shear  wave  velocity  from  the  shear-wave  value  of 
(»  800  ns)  measured  for  a  WC  test  flaw  and  the  diameter  estimates 
(1,160  urn  and  1,180  urn)  obtained  for  the  same  flaw  from  longitudinal-wave 
a  measurements.  This  resulted  in  the  desired  inverse-scattering  formula 

d  =  1.45  (10) 

where  d  is  in  micrometers  and  o ^  is  in  nanoseconds  with  an  effective  WC 
shear-wave  velocity  of  2.9  km/sec. 

2.6.3  Physical  Model  for  Sapphire  Inclusions 


Sapphire  Inclusions  are  nearly  impenetrable  due  to  the  magni 
tudes  of  their  LW  and  SW  velocities  (10  km/sec  and  6  km/sec,  respectively) 
relative  to  those  of  the  IN100  alloy  matrix.  Except  for  the  positive 
polarity  of  the  R  pulse,  the  impulse  response  of  a  sapphire  Inclusion  is 
expected  to  resemble  that  of  a  void  (see  the  lower  and  upper  diagrams  of 
Figures  7  and  8).  Equations  7  and  8  derived  for  voids  are  also  valid  for 
sapphire  Inclusions. 


3.  SOFTWARE  TRANSFER  AND  ADAPTATION 


3.1  General 


The  main  program  objective  was  to  independently  evaluate  the  Born 
inversion  software  developed  by  others  against  a  set  of  volumetric,  flaws 
under  blind  test  conditions.  To  accomplish  this  objective,  the  existing 
software  was  transferred  to  SwRI  and  made  operational  with  a  newly  assem¬ 
bled  hardware  system. 

3.2  Software  Transfer 


To  initiate  the  Born-inversion-software  adaptation  task,  a  subcon¬ 
tract  was  let  to  Ames  Laboratory  for  (1)  delivering  the  software  to  SwRI 
in  a  computer  language  suitable  for  the  inversion  of  deconvolved  and  prop¬ 
erly  time-shifted  waveforms  obtained  for  each  of  the  test  flaws,  and 
(2)  assisting  SwRI  in  making  the  Born  inversion  software  operational  with 
a  newly  assembled  data  acquisition  and  inversion  system. 

The  Ames  software  package  was  transferred  in  two  steps.  First,  the 
data-inversion  computer  program  (BIA)  for  estimating  the  flaw's  character¬ 
istic  function  and  effective  diameter  was  transferred  in  Fortran  language. 
In  the  Ames  data-processing  computer  program  (ZOTA),  the  origin  of  time  is 
chosen  at  the  maximal  cross-sectional  area.  The  software  transfer  was  com¬ 
pleted  in  November  1981. 

A  time-domain  deconvolution  algorithm  developed  at  the  Air  Force 
Institute  of  Technology  (4^)  was  also  transferred  to  SwRI  but  was  not  used 
extensively  after  its  equivalency  to  the  existent  Wiener-filter  type  decon¬ 
volution  algorithm  was  established  (23) . 

3.3  Data  Acquisition  and  Inversion  System 

The  hardware  comprising  the  Data  Acquisition  and  diversion  SYstem 
(DAISY)  assembled  for  specimen  testing  is  shown  in  Figure  9.  This  system 
is  composed  of  off-the-shelf  items  (ultrasonic  transducers,  pulsers, 
receivers,  etc.)  and  standard  Interfaces  and  connectors.  It  includes  a 
data  acquisition  stage,  a  data  storage  capability,  a  data  processing 
stage,  and  a  data  analysis  (Born  inversion)  stage.  These  are  described 
under  two  headings. 

3.3.1  Data  Acquisition  System 

The  ultrasonic  data  were  acquired  using  immersion  and  con¬ 
tact  transducers.  Because  of  the  limited  bandwidth  of  any  transducer, 
three  standard-frequency  wideband  transducers  (nominal  center  frequencies 
of  2.25,  5,  and  10  MHz)  were  made  available  so  that  at  least  one  of  the 
transducers  would  satisfy  the  0.5<ka<2.5  inequality  for  a  test  flaw  known 
to  be  in  the  250  to  1,300  um  size  range.  The  transducers  were  driven  by  a 
pulser-receiver  unit  made  by  Panametrics  (Model  5055PR)  that  enables  the 


transmission  of  short  (three  half-cycle)  ultrasonic  pulses  into  the  test 
specimen.  The  receiver  can  provide  up  to  60  dB  gain  and  68  dB  attenuation. 


A  digital  oscilloscope  made  by  Tektronics  (Model  7854)  was 
used  to  observe  the  received  waveforms  for  a  (satellite  pulse  separation) 
measurements.  The  waveform  analyzer  of  the  oscilloscope  digitizes  selected 
portions  (usually  2-ps  long  segments)  of  the  radio-frequency  (RF)  waveforms 
and  transfers  the  digitized  data  to  a  small  (desktop)  computer.  Each  wave¬ 
form  can  be  sampled  up  to  256  times  within  selectable  time  windows.  Many 
samples  were  acquired  and  then  averaged  to  minimize  the  background  inter¬ 
ference  caused  by  random  electronic  noise. 

3.3.2  Data  Inversion  System 


The  desktop  computer  for  the  DAISY  is  a  Hewlett-Packard  85 
and  the  main  computer  is  a  Data  General  Eclipse  C-330-  The  desktop  com¬ 
puter  offers  graphic  capability,  keyboard,  and  magnetic  tape  storage,  all 
in  a  single  portable  unit  located  in  the  ultrasonics  laboratory.  The  main 
functions  of  the  HP-85  are  to  store  the  digitized  waveforms  and  assist  in 
selecting  the  most  appropriate  transducer (s)  to  acquire  the  data  for  Born 
inversion.  The  software  package  for  the  HP-85  includes  programs  for  gene¬ 
rating  disk  files,  managing  the  data  storage  on  disk,  normalizing  the  flaw 
waveforms  with  the  reference  waveforms  (deconvolution),  and  displaying  the 
waveforms  and  the  complex  Fourier  spectra  (magnitude  and  phase  or  real  and 
imaginary  parts)  on  hard  copies. 

The  Born  inversion  software  from  Ames  laboratory  was  executed 
on  the  Eclipse  C-330  computer  along  with  two  deconvolution  programs  to 
ensure  timely  processing  and  Born  inversion  of  the  digitized  waveforms. 
Several  other  waveform  processing  programs  are  resident  in  this  computer, 
including  time  gating,  Fourier  transforming,  low-frequency  extrapolating, 
filtering,  and  time  shifting.  The  signals  in  the  digitized  flaw  waveforms 
(i.e,  the  main  signal  and  any  satellites  captured  in  a  record  length  of 
2  microseconds)  were  software  gated  and  Fourier  transformed.  The  flaw 
waveforms  were  then  deconvolved,  time  shifted,  and  put  into  the  Born  inver¬ 
sion  algorithm.  Hard  copies  of  the  Born  inversion  results  of  the  type  dis¬ 
played  in  Figure  6  were  obtained  for  each  test  flaw. 

3.4  Software  Adaptation 


The  SwRI  software-hardware  combination  was  tested  against  the  Ames 
combination  using  an  800  i.m  diameter  void  ("standard  flaw,"  2,4,23-25) . 
The  experimental  protocol  for  obtaining  the  Born  diameter  for  an  internal 
bulk  flaw  was  established  under  the  guidance  of  Ames  Laboratory  ( 1 , 26) . 
This  protocol  consists  of  the  steps  listed  in  Table  2.  Note  that  this 
protocol  did  not  concern  itself  with  a  method  for  selecting  the  trans¬ 
ducer^)  most  appropriate  for  a  flaw  of  unknown  type  and  size. 

The  "optimum"  transducer  for  acquiring  the  waveform  from  the  stan¬ 
dard  flaw  was  determined  by  the  Rockwell  transducer  selection  protocol  to 
be  the  5-MHz  transducer  0,2^4).  According  to  Equation  (1),  however,  the 
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"proper"  transducer  is  the  2.25  MHz  transducer.  Since  probably  neither  is 
the  "correct"  transducer  (27) ,  we  decided  to  apply  the  BIT  to  the  waveforms 
acquired  by  each  of  these  transducers.  The  Born  diameters  obtained  for  the 
standard  flaw  by  the  Ames  and  the  SwRI  hardware  systems  were  essentially 
the  same  (see  Table  3).  These  test  results  therefore  confirmed  that  the 
Born  inversion  software  was  successfully  adapted  by  SwRI. 


Table  2 


EXPERIMENTAL  PROTOCOL  FOR  THE  BORN  INVERSION  TECHNIQUE 


Description 


Digitize  and  store  the  flaw  waveform  y(t) 

Digitize  and  store  a  reference  waveform  w(t) 

Deconvolve  y(t)  with  w(t)  [see  Equation  (3)] 

Extrapolate  magnitude  spectrum  H(f)  to  low  frequencies 

Shift  impulse  response  h(t)  to  make  center  of  flaw 
coincide  with  origin  of  time 

Compute  characteristic  function  J(r)  [see  Equation  (4)] 
Estimate  flaw  size  [see  Equation  (5)] 

Repeat  steps  1  through  7  for  another  viewing  direction 


Table  3 

BORN  DIAMETER  ESTIMATES  OBTAINED  FOR  AN  800-um  DIAMETER  VOID 
AT  SOUTHWEST  RESEARCH  INSTITUTE  AND  AMES  LABORATORY 
WITH  TWO  TRANSDUCERS 


BORN  DIAMETER,  urn 


2.25  MHz 


5  MHz 


4.  PROCEDURE  DEVELOPMENT 


4.1  General 


The  primary  flaw  characterization  procedures  were  developed  based  on 
the  proven  capabilities  of  the  three  ultrasonic  techniques  and  the  three 
data  processing/inversion  algorithms  listed  in  Table  4.  These  techniques 
and  algorithms  are  described  in  Section  2.  The  spectrum  splicing  technique 
needs  supplementation  by  the  amplitude-  and  phase-comparison  techniques  and 
the  deconvolution  algorithm,  and  the  BIT  needs  supplementation  by  the  SST. 
Thus,  all  data  about  the  flaw  including  any  ji  priori  information  are  com¬ 
bined  to  produce  reliable  and  accurate  estimates  of  the  flaw  characteris¬ 
tics.  The  Born  inversion  procedures,  collected  under  three  headings,  spe¬ 
cifically  address  the  critical  transducer  selection  problem  for  a  flaw  of 
unknown  type  and  size. 

4.2  Flaw  Identification  Procedures 


Preliminary  examination  of  the  "raw"  composite  waveforms  obtained 
with  the  high-frequency  (5  and  10  MHz)  transducers  for  the  spatial  pat¬ 
terns  and  amplitude  and  phase  distributions  of  the  constituent  pulses  can 
reveal  the  likely  flaw  composition  and  shape.  If  the  "locomotive”  of  the 
pulse  train  (i.e.,  the  R  pulse  originating  from  the  front  surface  of  the 
flaw)  has  a  positive  polarity  and  If  the  pulses  that  follow  closely  are  of 
comparable  amplitude  and  like  polarity,  then  the  flaw  is  likely  to  be  a 
weak  inclusion  (e.g.,  tungsten  carbide).  If  the  R  pulse  with  a  positive 
polarity  is  not  followed  by  a  discernible  pulse  of  like  polarity,  then  the 
flaw  is  likely  to  be  a  strong  inclusion  (e.g.,  sapphire).  Finally,  if  the 
R  pulse  has  a  negative  polarity,  then  the  flaw  is  likely  to  be  a  void. 


Table  4 

SOURCES  AND  EXPECTED  RESULTS/OUTPUTS  OF  TECHNIQUES/ 
ALGORITHMS  INCORPORATED  INTO  THE  PRIMARY 
FLAW  CHARACTERIZATION  PROCEDURES 


Technique /Algorithm 

Result/Output 

Source 

Amplitude  comparison 

Phase  comparison 

Likely  flaw  composition 
and  shape 

Unknown 

Deconvolution 

Magnitude  spectra 

SwRI , 

Air  Force 

Spectrum  splicing 

Appropriate  transducer (s ) 

Ames , 
SwRI 

Rockwell , 

Zero-of-t ime  ^ ) 

Shifted  impulse  response 

Ames 

Born  inversion^) 

Flaw  size  estimate 

Ames 

^■These  algorithms  are  the 

essence  of  the  Born  inversion 

technique , 

4.3  Transducer  Selection  Procedures 


Qualitative  examination  of  the  deconvolved  flaw  spectrum  H(f) 
obtained  with  a  LW  transducer  for  its  extent  of  modulation  can  reveal  the 
appropriateness  of  that  transducer  for  inclusion  into  the  Born  inversion 
protocol  (Table  2).  Four  transducers  were  used  to  obtain  a  set  of  decon¬ 
volved  spectra  for  each  of  three  spherical  voids  of  known  size  in  an  effort 
to  test  the  ability  of  the  SST  to  quickly  rule  a  transducer  "appropriate" 
or  "inappropriate"  for  that  flaw.  The  spectra  obtained  for  a  550-i.m  diam¬ 
eter  void  in  an  IN100  specimen  are  shown  in  Figure  10  against  the  back¬ 
grounds  of  the  individual  transducer  spectra  (dashed  curves).  The  orders 
of  the  spectral  (resonance)  peaks,  defined  in  Figure  3,  can  be  readily 
identified.  The  center  frequencies  of  the  10-MHz  probes  are  distant  from 
the  flaw's  fundamental  (n=l)  resonance  frequency.  These  probes  are  then 
inappropriate  for  acquiring  Born  inversion  data  for  this  flaw.  Note  that 
the  center  frequency  of  the  5-MHz  probe  lies  to  the  right  of  the  first 
spectral  peak  and  that  of  the  2-MHz  probe  lies  to  the  left  of  the  first 
spectral  peak.  Both  of  these  probes  are,  therefore,  appropriate  for 
acquiring  Born  inversion  data  for  this  flaw.  Note  that  the  "proper"  trans¬ 
ducer  with  a  center  frequency  equal  to  the  average  spectrum  periodicity 
(3.4  MHz)  was  not  available. 

The  magnitude  spectra  obtained  for  the  800-um  diameter  (standard) 
void  by  sequentially  acquiring  waveforms  with  transducers  having  succes¬ 
sively  lower  center  frequencies  are  shown  in  Figure  11.  The  2. 5-MHz  trans¬ 
ducer  whose  bandwidth  is  almost  perfectly  centered  on  the  first  peak  of 
H(f)  is  the  most  appropriate  one  for  acquiring  the  backscattered  waveforms 
for  Born  inversion. 

As  expected,  the  smallest  spectrum  periodicities  (p  values)  were  mea¬ 
sured  for  the  largest  void  (see  Figure  12).  While  the  high-frequency 
transducers  served  quite  well  in  determining  the  Born  frequency  for  this 
strong  flaw  (F  =  2  MHz),  none  of  these  transducers  is  appropriate  for 
acquiring  the  waveforms  for  Born  inversion. 

The  spectral  peaks  can  be  more  reliably  ordered  when  the  individual 
spectra  obtained  for  the  same  flaw  are  spliced  together  (see  Figure  13). 
Note  that  for  only  one  of  the  three  flaws  did  the  required  frequency  F  (=p) 
coincide  with  the  effective  center  frequency  of  one  of  the  three  standard- 
frequency  transducers  (i.e.,  for  the  largest  flaw). 

The  purpose  of  the  transducer  selection  protocol  is  to  rule  at  least 
one  of  the  three  standard-frequency  transducers  inappropriate  for  acquir¬ 
ing  the  Born  inversion  data.  The  gates  (time  windows)  for  the  SST  wave¬ 
forms  must  include  all  trailers  of  the  R  pulse  recognized  during  the  quali¬ 
tative  flaw-identification  process.  The  protocol  for  selecting  the  most 
appropriate  Born  transducer(s)  consists  of  the  steps  listed  in  Table  5. 

4.4  Flaw  Sizing  Procedures 


A  subcontract  was  let  to  Rockwell  Science  Center  to  assure  that  the 
Ames  sizing  procedures  featuring  the  BIT  are  basically  the  same  as  the 
Rockwell  sizing  procedures  (i.e.,  they  yield  essentially  the  same  size 
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Figure  12.  Magnitude  Spectra  Obtained  for  a  1,050-um  Diameter  IN100 
Void  With  Four  Longitudinal-Wave  Transducers.  Note  that 
the  fundamental  resonance  frequency  of  this  flaw  is  2.0  MHz. 
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Figure  13.  Application  of  the  Spectrum  Splicing  Technique 


Table  5 


EXPERIMENTAL  PROTOCOL  FOR  THE  SPECTRUM  SPLICING  TECHNIQUE 

Step _ Description _ 

1  Set  the  gate  for  the  flaw  waveform  y(t)  obtained  with  the 
10-MHz  transducer 

2  Digitize  and  store  y(t) 

3  Digitize  and  store  a  gated  reference  waveform  w(t) 

4  Deconvolve  y(t)  with  w(t)  [see  Equation  (3)] 

5  Repeat  steps  1  through  4  with  the  5-MHz  transducer 

6  Repeat  steps  1  through  4  with  the  2.25-MHz  transducer 

7  Splice  the  individual  magnitude  spectra  H(f)  together  (see 
Figure  13) 

8  Determine  average  spectrum  periodicity  p 

9  Calculate  required  frequency  F  [see  Equation  (2)] 

10  Eliminate  the  inappropriate  transducer (s)  from  the  Born 
inversion  protocol  (Table  2) 


estimates  for  the  same  set  of  stored  waveforms).  The  standard  flaw  was 
the  basis  of  comparison.  Figures  14  through  17  contain  the  Born  plots 
H(f),  h(t),  A(f),  and  J(r)  obtained  by  applying  the  Ames  sizing  procedures 
( 1 , 13,26)  to  the  waveforms  acquired  at  SwRI  with  transducers  having  suc¬ 
cessively  lower  center  frequencies.  The  amplitude  spectrum  A(f)  and  the 
characteristic  function  J(r)  form  a  transform  pair  [see  Equation  (4)]. 

The  deconvolved  flaw  spectrum  H(f)  and  the  impulse  response  h(t)  are  also 
known  to  form  a  transform  pair.  It  can  be  shown  (13)  that  the  periodicity 
of  the  flaw's  magnitude  spectrum  is  governed  by  the  extent  of  its  impulse 
response  (i.e.,  p  =  l/o).  A  similar  reciprocal  relationship  appears  to 
exist  between  the  width  of  the  characteristic  function  and  the  mean  fre¬ 
quency  of  the  amplitude  spectrum  (i.e.,  a*  =  1/k*).  The  mean  frequency 
is,  in  turn,  governed  by  the  transducer  center  frequency.  Other  param¬ 
eters  being  the  same,  the  larger  the  transducer  center  frequency,  the 
smaller  the  Born  diameter  (see  Table  6). 

Figures  18  through  21  contain  the  Born  plots  H(f),  D(x),  and  J(t) 
obtained  by  applying  the  Rockwell  sizing  procedures  (2^)  to  the  waveforms 
acquired  either  at  SwRI  or  Rockwell  with  transducers  having  successively 
lower  center  frequencies.  The  main  difference  between  the  Rockwell  and 
Ames  Born  inversion  procedures  is  in  the  way  the  center  of  the  flaw  is 
located.  The  Ames  ZOTA  is  described  in  connection  with  Figures  5  and  6. 
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Figure  15.  Transform  Pairs  Obtained  for  an  800-»m  Diameter  Void  With 
a  4.2  MHz  Longitudinal-Wave  Transducer  Using  the  Ames 
Sizing  Procedures.  Note  that  po  =  0.96  and  k*a*  *  0.91. 


0 


10  MHz 


Figure  16.  Transform  Pairs  Obtained  for  an  800-um  Diameter  Void  With 
a  3.6  MHz  Longitudinal -Wave  Transducer  Using  the  Ames 
Sizing  Procedures.  Note  that  po  =  1.13  and  k*a*  =*  0.90. 


MAGNITUDE  SPECTRUM 


Figure  18.  Born  Diameters  Obtained  for  an  800-um  Diameter  Void  With 
a  6.2  MHz  Longitudinal-Wave  Transducer  Using  the  Rockwell 
Sizing  Procedures.  Note  that  kna  =  1.20. 


CHARACTERISTIC  FUNCTION 
FOR  t  =  96  ns 


Figure  19.  Born  Diameters  Obtained  for  an  800-;. m  Diameter  Void  With 
a  4.3  MHz  Longitudinal-Wave  Transducer  Using  the  Rockwell 
Test  Bed  System  (Reference  2)  and  Sizing  Procedures. 

Note  k0aQ  =  1.31. 
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VARIATION  OF  BORN  DIAMETER 
WITH  TIME  SHIFT  t 


dQ  =  480  ym 


CHARACTERISTIC  FUNCTION 
FOR  t  =  10  ns 


Born  Diameters  Obtained  for  an  800-um  Diameter  Void  With 
a  4.2  MHz  Longitudinal-Wave  Transducer  Using  the  Rockwell 
Sizing  Procedures.  Note  that  kQa0  =  1.14. 


Table  6 


BORN  DIAMETER  ESTIMATES  OBTAINED  WITH  FIVE  TRANSDUCERS 
FOR  AN  800-um  DIAMETER  VOID  AT  SOUTHWEST.  RESEARCH 
INSTITUTE  AND  ROCKWELL  SCIENCE  CENTER 


Transducer  Center 
Frequency,  MHz 
Nominal  Effective 


SwRI  (oe,  ns)  (4) 


Born  Diameter, 
Rockwell 


d 


o 


10 

6.2 

290 

(370) 

480 

280 

5 

4.2 

490 

(370) 

740 

480 

5 

4.3(2) 

— 

( - ) 

760 

510 

5 

3.6 

650 

(370) 

— 

— 

2.25 

2.5(3) 

790 

(370) 

1,250 

800 

Size 

estimates  obtained 

with  the 

most  appropriate 

(SST-approved) 

trans- 

ducers  are  underlined. 

2  , 

"Optimum"  transducer  per  Rockwell  procedures  (Reference  2). 
^"Proper"  transducer  per  SwRI  procedures  [Equation  (1)]. 

4cre=  1/p  (Figure  11). 


The  Rockwell  ZOTA  has  no  physical  basis.  An  empirical  technique  has  been 
found  (2)  which  works  when  the  transducer's  bandwidth  is  well  centered  on 
the  first  peak  of  the  magnitude  spectrum  (see  upper  diagram  in  Figure  21). 
This  technique  requires  that  the  characteristic  function  J(t)  have  maximum 
flatness  at  t=0. 

Two  steps  are  involved  in  obtaining  the  required  time  shift  for  the 
impulse  response.  First,  the  Born  diameter  is  determined  for  a  wide  range 
of  time  shifts  (see  middle  diagram  of  Figure  19).  It  is  known  that  the 
center  of  the  flaw  is  located  near  the  second  maximum  of  D(t)  and  dis¬ 
placed  to  the  left  of  it.  Next,  a  series  of  characteristic  functions  is 
generated  for  time  shifts  ranging  from  t=0  to  100  ns.  The  characteristic 
function  was  maximally  flat  for  t  =96  ns  and  it  resulted  in  a  Born  diameter 
of  760  ,  m  (see  lower  diagram  in  Figure  19).  However,  in  another  case  when 
the  transducer  bandwidth  was  also  well  centered  on  the  first  peak  of  the 
magnitude  spectrum  (see  upper  diagram  of  Figure  21),  the  Born  diameter 
corresponding  to  the  maximally  flat  characteristic  function  was  1,250  ,,m. 

The  Rockwell  column  in  Table  6  contains  the  Born  diameters  obtained 
by  maximizing  the  flatness  of  the  characteristic  function  at  t=0.  Rather 
than  these  estimates,  the  Born  diameters  denoted  by  d0  in  Figures  18 
through  21  and  collected  in  the  last  column  of  Table  6  appear  to  correspond 
to  those  obtained  with  the  Ames  procedures  and  included  in  the  SwRI  column 
of  Table  6.  This  established  the  equivalency  of  the  Ames  and  Rockwell  Born 
inversion  procedures. 
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5.  TEST  SPECIMENS 


5.1  General 


In  order  to  qualify  the  developed  flaw  characterization  procedures 
it  is  necessary  to  have  test  specimens  that  contain  a  statistically  large 
number  of  internal  bulk  flaws  with  characteristics  (other  than  locations) 
unknown  to  the  ultrasonic  examiners.  A  subcontract  was  let  to  Pratt  & 
Whitney  Aircraft  Group  (PWA)  for  the  fabrication  of  titanium  alloy  and 
powder  metallurgy  nickel-based  superalloy  test  specimens  containing  40 
simulated  flaws  and  the  collection  of  up  to  10  real-world  flaws.  Flaw 
locations  and  characteristics  (composition,  shape,  size,  and  orientation) 
were  documented  by  PWA  during  specimen  design.  However,  only  the  flaw 
locations,  the  250  to  1,300  um  size  range,  and  the  list  of  possible  flaw 
compositions  were  submitted  to  SwRI  with  the  specimens  to  preserve  the 
credibility  of  the  test  results.  This  section  describes  the  material  and 
flaw  types,  the  test  specimens,  and  the  procedures  employed  to  produce 
these  specimens. 

5.2  Material  and  Flaw  Types 


Titanium  alloys  and  nickel-based  superalloys  are  the  two  basic  mate¬ 
rial  types  used  to  make  rotor  components  for  gas  turbine  engines.  In  each 
alloy,  flaws  of  different  composition  have  been  found  according  to  the 
differences  in  the  chemical  activity  of  the  alloy  base  elements  and  in  the 
manufacturing  processes  used  to  form  these  materials  into  engine  disks. 

Three  types  of  manufacturing  flaws  have  been  identified  in  titanium 
alloys — porosity,  chemical  segregation,  and  tungsten  carbide  (WC)  inclu¬ 
sion.  Voids  and  WC  inclusions  are  readily  detectable  with  ultrasonic 
waves.  Manufacturing  flaws  are  found  more  frequently  in  nickel  alloys 
than  in  titanium  alloys  because  melt  practice  is  not  yet  sufficiently 
developed  to  prevent  the  introduction  of  voids  and  hard  oxide  inclusions 
(e.g.,  sapphire). 

Disk  material  specimens  containing  46  separately  identifiable  voids 
and  inclusions  were  provided  by  PWA  of  West  Palm  Beach,  Florida  (28) .  The 
shapes  of  the  40  electro-discharge  machined  (EDM)  flaws  were  limited  to 
spheres.  Flaw  sizes  were  to  be  in  the  250  to  1,300  ..m  range.  The  distri¬ 
bution  of  the  disk  materials,  flaw  types,  and  number  of  flaws  per  specimen 
is  shown  in  Table  7.  Inclusions  in  the  diffusion-bonded  titanium  specimen 
were  to  be  tungsten  carbide  inclusions  and  those  in  the  diffusion-bonded 
IN100  specimens  were  to  be  sapphire  (AI2O3)  inclusions. 

Manufacturing  processes  (machining,  bonding,  heat  treating,  and  sur¬ 
face  finishing)  on  the  test  specimens  were  the  same  as  on  turbine  engine 
components.  Special  attempts  were  made  to  allow  for  a  clear  line  of  sight 
to  the  test  flaws  from  all  sides  of  the  diffusion-bonded  specimens. 


Table  7 

TEST  SPECIMENS  PROVIDED  BY  PRATT  &  WHITNEY  AIRCRAFT  GROUP 


Specimen  No. 

Material 

Flaw  Type 

No.  of  Flaws  (ID  Nos.) 

FML  97413 

Ti  6-4 

Void 

5  (01-05) 

FML  97414 

IN100 

Void 

5  (06-10) 

FML  97416 

Ti  6-4 

Inclusion 

5  (11-15) 

FML  97474-1 

IN100 

Inclusion 

3  (16-18) 

-2 

IN100 

Inclusion 

3  (19-21) 

-3 

IN100 

Inclusion 

3  (22-24) 

-4 

IN100 

Inclusion 

3  (25-27) 

-5 

IN100 

Inclusion 

3  (28-30) 

-6 

IN100 

Inclusion 

3  (31-33) 

-7 

INI  00 

Inclusion 

4  (34-37) 

-8 

IN  100 

Inclusion 

3  (38-40) 

FML  97600-1 

IN100 

Unknown 

3  (41-43) 

-2 

3  (44-46) 

Total 


5.3  Voids  In  Titanium  Specimen  No.  FML  97413 

This  specimen  consists  of  three  diffusion-bonded  plates  of  varying 
thickness.  Five  spherical  voids  with  diameters  ranging  from  320  to 
1,230  um  are  located  in  the  two  bondplanes.*  The  manufacturing  processes 
for  this  specimen  are  listed  in  the  second  column  of  Table  8. 

5.4  Inclusions  in  Titanium  Specimen  No.  FML  97416 


This  specimen  is  configured  the  same  way  as  specimen  No.  FML  97413. 
Five  spherical  WC  inclusions  with  diameters  ranging  from  480  to  1,210  wm 
are  located  in  the  two  bondplanes.  The  manufacturing  processes  were  dif¬ 
ferent  for  this  titanium  specimen  (see  third  column  in  Table  8).  A  hot 
isostatic  press  (HIP)  process  was  used  to  close  the  holes  around  the  inclu 
sions  after  the  plates  were  bonded  together.  This  was  done  to  ensure  con¬ 
tinuous  contact  between  the  host  material  and  the  WC  spheres. 

5 . 5  Voids  in  IN100  Specimen  No.  FML  97414 

This  specimen  consists  of  three  super  powder  metallurgy  IN100  plates 
diffusion-bonded  with  nickel  foil.  Five  spherical  voids  with  diameters 

*The  ultrasonic  examiners  did  not  know  that  the  voids  were  spherical. 
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Table  8 


MANUFACTURING  PROCESS  PLANS  FOR 

Sequence  Void  Specimen  No.  FML  97413 

1  Cut  raw  material  from  stock 

2  Rough  cut  and  drill  two 
position  holes 

3  Grind  and  lap  bonding 
surfaces 

4  Mark  plates 

5  EDM  voids  (5) 

6  Deburr  holes  by  light  lap 

7  Fabricate  two  position  pins 

8  Assemble  and  diffusion-bond 
plates 

9  Heat  treat  to  Spec.  PWA  1215 

10  Final  machine 

11  Ultrasonic  examination 

12  Ship  to  SwRI 

13 


TITANIUM  6-4  TEST  SPECIMENS 

Inclusion  Specimen  No.  FML  97416 

Cut  raw  material  from  stock 

Rough  cut  and  drill  two  position 
holes 

Grind  and  lap  bonding  surfaces 
Mark  plates 

EDM  holes  for  inclusions  (5) 

Deburr  holes  by  light  lap 
Fabricate  two  position  pins 
Assemble  and  diffusion-bond  plates 

Hot  isostatic  press  (HIP) 

Heat  treat  to  Spec.  PWA  1215 
Final  machine 
Ultrasonic  examination 
Ship  to  SwRI 


ranging  from  320  to  1,230  urn  are  located  in  the  two  bondplanes.  The  manu¬ 
facturing  processes  for  this  specimen  are  listed  in  the  second  column  of 
Table  9. 


5.6  Inclusions  in  IN100  Specimen  Nos.  FML  97474-1  through  -8 

Eight  IN100  specimens  containing  a  total  of  25  spherical  sapphire 
inclusions  were  fabricated.  Each  test  specimen  is  made  up  of  two  diffusion- 
bonded  plates  of  equal  thickness.  The  inclusions  with  diameters  ranging 
from  410  to  990  urn  are  located  in  the  bondplane.  The  manufacturing  pro¬ 
cesses  for  these  specimens  are  listed  in  the  third  column  of  Table  9- 

5. 7  Real  Flaws  in  IN100  Specimen  Nos.  FML  97600-1  and  -2 

Two  IN100  specimens  were  provided  with  six  flaws  of  unknown  composi¬ 
tion,  shape,  size,  and  orientation  found  during  production  material  pro¬ 
cessing  in  an  extruded  billet.  The  billet  was  forged  into  a  pancake  from 
which  the  specimens  were  cut. 
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Table  9 


MANUFACTURING  PROCESS  PLANS  FOR  IN100  TEST  SPECIMENS 


Sequence  Void  Specimen  No.  FML  97414 

1  Cut  raw  material  from  stock 

2  Rough  cut  and  drill  two 
position  holes 

3  Grind  and  lap  bonding 
surfaces 

4  Mark  plates 

5  EDM  Voids  (5) 

6  Deburr  holes  by  light  lap 

7  Fabricate  two  position  pins 

8  Assemble  and  diffusion-bond 
plates 

9  Heat  treat  to  modified  Spec. 
PWA  1074 

L0  Final  machine 

11  Ultrasonic  examination 

12  Ship  to  SwRI 

13 


Inclusion  Specimen  No.  FML  97474 
Cut  raw  material  from  stock 
Rough  cut  plates 

Grind  and  lap  bonding  surfaces 

EDM  holes  for  inclusions  (3-5) 
Deburr  holes  by  light  lap 
Mark  plates 

Assemble  and  diffusion-bond  plates 
HIP 

Heat  treat  to  modified  Spec. 

PWA  1074 

Final  machine 

Mark  specimen 

Ultrasonic  examination 

Ship  to  SwRI 
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6.  PROCEDURE  QUALIFICATION 


6.1  General 


The  simulated  and  real  flaws  in  the  test  specimens  were  used  solely 
for  verification  of  the  applicability  of  the  selected  flaw  characterization 
technique  combinations  described  in  Section  2.  In  order  to  reconstruct 
the  flaw  from  the  results  of  the  one-dimensional  BIT  it  was  necessary  to 
acquire  many  flaw  waveforms  from  regularly  placed  transducer  positions. 
While  the  flaws  were  accessible  from  all  six  sides  of  a  test  specimen,  no 
measurements  were  taken  from  the  four  sides  intersecting  the  bondplanes 
associated  with  the  specimen  manufacturing  process. 

The  likely  compositions,  shapes,  and  sizes  for  the  42  detectable  test 
flaws  are  presented  in  four  parts.  The  first  two  parts  (Subsections  6.2 
and  6.3)  contain  the  results  obtained  by  the  serial  applications  of  the 
supplemented  SST  and  the  BIT  [i.e.,  the  techniques  incorporated  into  the 
primary  flaw  characterization  procedures  (see  Table  4)].  The  SPOT  flaw 
size  estimates  are  given  in  Subsection  6.4.  Finally,  the  results  obtained 
for  the  32  detectable  inclusions  by  an  automatable  version  of  the  SPOT-SST 
combination  that  waives  the  need  for  a  highly  trained  operator  are  pre¬ 
sented  in  Subsection  6.5. 

6.2  Born  Frequency  Estimates 

The  transducer  selection  results  obtained  by  the  application  of  the 
SST  to  the  individual  magnitude  spectra  of  each  test  flaw  are  presented 
under  four  headings  (i.e.,  voids,  tungsten  carbide  inclusions,  sapphire 
inclusions,  and  real  flaws). 

6.2.1  Voids  in  Titanium  and  IN100  Specimens 

The  amplitude-  and  phase-comparison  techniques  were  applied  first  to 
the  flaws  in  titanium  specimen  No.  FML  97413  and  IN100  specimen  No.  FML 
97414.  As  the  10-MHz  (or  5-MHz)  immersion  LW  transducer  was  moved  around 
a  given  flaw,  the  received  waveforms  were  noted  to  contain  essentially  the 
same  key  signal  parameters  (pulse  polarities,  pulse  train  durations,  and 
echo  dynamic  patterns).  Based  on  these  observations,  the  flaws  in  speci¬ 
men  Nos.  FML  97413  and  FML  97414  were  identified  to  be  spherical  voids.* 

The  individual  magnitude  spectra  H(f)  obtained  for  one  of 
the  voids  (flaw  02)  with  two  10-MHz,  one  5-MHz,  and  one  2-MHz  transducers 
are  displayed  against  the  backgrounds  of  the  appropriate  transducer  spec¬ 
tra  W(f)  in  Figure  22.**  The  first  five  external  resonances  (spectral 
peaks  1  through  5)  of  this  flaw  were  readily  identified  in  each  of  the 


*These  preliminary  conclusions  were  later  substantiated  by  SW  measurements 
(see  Subsection  6.4). 

**Appendix  A  contains  a  full  set  of  magnitude  spectra. 


individual  plots.  However,  the  first  two  spectral  peaks  of  flaw  01  were 
not  identified  until  the  individual  spectra  obtained  with  the  two  10-MHz 
transducers  were  spliced  together  (see  Figure  23).  The  improved  visibil¬ 
ity  of  the  spectral  peaks  is  the  merit  of  the  SST.  The  flaw  profiles 
(spliced  frequency-domain  responses)  of  Figure  23  indicate  that  the  2-MHz 
transducer  is  the  proper  one  for  flaws  02  and  05,  while  the  5-MHz  and  the 
10-MHz  transducers  should  be  used  in  acquiring  the  Born  data  on  the  other 
three  flaws.  The  profiles  obtained  for  flaws  06  through  10  are  reproduced 
in  Figure  24.  The  effective  transducer  center  frequencies  (F)  required  by 
the  ten  voids  are  listed  in  the  last  column  of  Table  10. 

6.2.2  Tungsten  Carbide  Inclusions  in  Titanium  Specimen 

The  application  of  the  amplitude-  and  phase-comparison  tech¬ 
niques  to  the  flaws  in  titanium  specimen  No.  FML  97416  revealed  that  these 
flaws  are  spherical  tungsten  carbide  inclusions.*  Figure  25  shows  the 
individual  magnitude  spectra  obtained  for  one  of  the  WC  inclusions  (flaw 
11)  with  four  immersion  LW  transducers.**  The  profiles  obtained  for  the 
five  WC  inclusions  are  shown  in  Figure  26.  The  effective  transducer  center 
frequencies  required  by  these  inclusions  are  obtained  by  multiplying  the 


Table  10 

BORN  FREQUENCY  ESTIMATES  OBTAINED  FOR  THE  VOIDS  IN  TITANIUM 
SPECIMEN  NO.  FML  97413  AND  IN100  SPECIMEN  NO.  FML  97414  BY 
THE  SPECTRUM  SPLICING  TECHNIQUE  (SST) 


Flaw 
ID  No. 


_ Transducer  Center  Frequency _ 

Low  Med .  Low  Med.  High  High 


Average  Spectral 
Period,  MHz^1) 


01 

—  (2) 

—  (2) 

6.5 

5.6 

6.0 

02 

2.0 

2.0 

2.0 

2.1 

2.0 

—  (2) 

3.2 

2.7 

2.5 

2.8 

04 

—(2) 

2.8 

2.5 

2.5 

2.6 

05 

2.0 

1.7 

2.0 

2.0 

1.9 

06 

2.7 

3.0 

2.7 

2.8 

1.2 

1.4 

—  (3) 

1.3 

08 

—(2) 

5.4 

6.4 

5.9 

09 

—  (2) 

2.1 

2.3 

2.2 

10 

2.5 

2.1 

2.4 

2.3 

lThe  Born  frequency  Cor  a  void  is  given  by  the  average  periodicity  of  its 
magnitude  spectra  (spectral  period,  F*=p). 

^Magnitude  spectrum  peaked  at  the  transducer  center  frequency  (feedthrough) 
■^Unmodulated  magnitude  spectrum. 


*These  preliminary  conclusions  were  borne  out  by  SW  measurements  (see 
Subsection  6.4). 

**Appendix  A  contains  a  full  set  of  magnitude  spectra. 
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•  Figure  23.  Application  of  the  Spectrum  Splicing  Technique  (SST) 

^  to  the  Voids  in  Titanium  Specimen  No.  FML  97413 

[■  (Flaws  01  Through  05) 
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de  Spectra  Obtained  for  Titanium  Inclusion 
With  Four  Longitudinal-Wave  Transducers 


average  spectral  periods  listed  In  the  last  column  of  Table  11  by  0.7 
[see  Equation  (2)]. 

6.2.3  Sapphire  Inclusions  In  IN100  Specimens 

Due  mainly  to  the  low  impedance  mismatch  between  the  flaws 
and  host  material  in  the  IN100  specimen  Nos.  FML  97474-1  through  -8,  four 
out  of  the  total  of  25  flaws  were  not  detected  at  all,  and  only  five  flaws 
were  detected  with  more  than  one  immersion  LW  transducer.  The  low  signal 
amplitudes,  the  positive  polarities  of  the  specularly  reflected  signals, 
the  absence  of  internal  echoes,  and  the  invariance  of  the  main  features  of 
the  flaw  waveform  with  specimen  rotation  have  led  us  to  conclude  that  these 
flaws  are  spherical  sapphire  inclusions.*  Figure  27  contains  the  profiles 
obtained  for  the  five  sapphire  inclusions  that  were  detected  more  than 
once.  The  effective  transducer  center  frequencies  required  by  these  inclu¬ 
sions  are  listed  in  the  last  column  of  Table  12. 

6.2.4  Real  Flaws  in  IN100  Specimens 

The  SST  could  not  be  applied  to  the  manufacturing  flaws  in 
INL00  specimen  Nos.  FML  97600-1  and  -2  because  only  one  immersion  LW  trans¬ 
ducer  was  effective  in  detecting  these  flaws.  The  magnitude  spectra 
obtained  for  flaws  41  through  46  with  this  transducer  are  collected  in 
Figure  A-19. 

6.3  Born  Diameter  Estimates 


Originally,  the  BIT  was  to  be  applied  only  to  the  waveforms  obtained 
with  a  transducer  appropriate  for  the  flaw  (see  the  last  columns  in  Tables 

Table  11 

BORN  FREQUENCY  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN 
TITANIUM  SPECIMEN  NO.  FML  97416  BY  THE  SST 


Flaw  _ Transducer  Center  Frequency _ 

ID  No .  Low  Med.  Low  Med.  High  High 


Average  Spectral 
Period,  Mliz^1) 


11 

—  (2) 

3.3 

3.7 

3.5 

3.5 

12 

2.7 

2-7 

2.6 

2.5 

2.6 

13 

2.6 

3.0 

2.7 

2.5 

2.7 

14 

— 

3.2 

3.9 

3.5 

3.5 

15 

— 

6.0 

7.0 

5.7 

6.2 

^■The  Born  frequency  for  a  weak  inclusion  (Born  flaw)  is  given  by  70  percent 
of  the  average  spectral  period  (F=0.7p). 

''Magnitude  spectrum  peaked  at  the  transducer  center  frequency  (feedthrough). 


*The  same  conclusions  were  reached  later  by  SW  measurements  (see  Subsec¬ 
tion  6.4). 
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Table  12 


BORN  FREQUENCY  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN 
IN  100  SPECIMEN  NOS.  FML  97474-3  AND  -4  BY  THE  SST 


Flaw 

Transducer 

Center 

Frequency 

Average  Spec 

ID  No. 

Low 

Medium 

High 

Period,  Mllz 

22 

2 . 2 

2.2 

2.3 

2.2 

23 

2.4 

2.3 

2.1 

2.3 

24 

2.2 

2.2 

2.2 

2.2 

2  5 

—  (2) 

3.3 

3.2 

3.3 

27 

— 

2.5 

2.5 

2.5 

^The  Born  frequency  for  a  strong  inclusion  (Franz  flaw)  is  given  by  the 
average  spectral  period  (F=p). 

“^Magnitude  spectrum  peaked  at  the  effective  transducer  center  frequency 
(feedthrough) . 


10  through  12  for  the  desired  effective  transducer  center  frequencies). 
However,  for  most  of  the  flaws  a  transducer  of  the  required  effective 
center  frequency  was  not  available  or  efficient  .  It  was  also  recognized 
that  a  large  inadvertent  error  in  the  SST  results  would  lead  to  a  corre¬ 
sponding  error  in  the  Born  diameter  by  forcing  the  choice  of  an  incorrect 
transducer  ( 2_7 ) .  For  these  reasons,  the  data  acquired  for  a  flaw  with  a 
seemingly  inappropriate  transducer  were  also  inverted  while  still  working 
under  blind  test  conditions.  These  data  are  included  to  illustrate  the 
criticality  of  the  transducer  selection  process.  To  set  the  size  esti¬ 
mates  obtained  with  inappropriate  transducers  apart  from  those  obtained 
with  the  SST-approved  transducers,  in  the  tables  that  follow  the  latter 
results  are  underlined-  The  results  are  presented  separately  for  voids, 
WC  inclusions,  sapphire  inclusions,  and  real  flaws. 

6.3.1  Voids  in  Titanium  and  1N100  Specimens 


The  desired  effective  transducer  center  frequencies  for  the 
ten  voids  are  reproduced  from  Table  10  for  ready  reference  in  the  second 
column  of  Table  13.  Listed  also  in  Table  13  are  the  nominal  flaw  diameters 
published  by  PWA  (2_8)  for  comparison  with  the  Born  diameters  obtained  with 
three  standard-frequency  transducers.  Note  that  the  requirement  made  on 
i he  effective  transducer  center  frequency  by  flaw  01  (i.e.,  6  MHz)  could 
not  be  met  by  any  of  the  available  2-MHz,  4-MHz,  ind  8-MHz  transducers. 

For  this  reason,  two  Born  diameters  (240  ,  m  and  470  .  m)  are  underlined  for 
this  flaw.  The  nominal  flaw  diameter  (320  m)  is  properly  bracketed  by 
the  underlined  estimates  because  the  SST  estimate  for  F  (6  MHz)  was  close 
to  its  expected  value  of  7  MHz  [=2,200/320  per  Equation  (1)]. 


Table  13 


BORN  DIAMETER  ESTIMATES  OBTAINED  FOR  THE  VOIDS  IN  TITANIUM 
SPECIMEN  NO.  FML  97413  AND  IN100  SPECIMEN  NO.  FML  97414 
WITH  THREE  LONGITUDINAL-WAVE  TRANSDUCERS 


Flaw 

ID  No. 

Required 
Freq .  ,  MHz 

Born 

2  MHz  ~ 

Diameter . 
4  MHz 

jj  m(  1 ) 

IT  MHz 

Nominal  Diam 

01 

6.0 

1,000 

470 

240 

320 

02 

2.0 

1,030 

450 

250 

1,000 

03 

2.8 

990 

510 

250 

600 

04 

2.6 

960 

530 

240 

780 

05 

1.9 

1,000 

560 

230 

1,230 

06 

2.8 

800 

540 

240 

600 

07 

1.3 

970 

530 

210 

1,230 

08 

5.9 

1,050 

440 

230 

320 

09 

2.2 

910 

480 

240 

780 

10 

2.3 

910 

440 

220 

1,000 

^Size  estimates  obtained  with  the  most  appropriate  transducers  are 
underlined . 


Appendix  B  contains  representative  sets  of  amplitude  spectra 
(inputs  to  the  BIA)  and  characteristic  functions  (outputs)  from  which  the 
Born  diameters  D  (and  d*)  are  derived  in  accordance  with  Equation  (5). 

Note  that  the  narrower  the  amplitude  spectrum  (i.e.,  the  smaller  the  mean 
frequency  f*),  the  broader  the  characteristic  function  (i.e.,  the  larger 
the  Born  diameter  d*)  and  vice  versa  (see  Figure  B-l  for  the  definitions 
of  f*  and  d*).  Note  that  the  average  value  of  the  k*a*  product  is  0.91. 

Flaw  02  is  of  special  interest  in  that  a  transducer  whose 
effective  center  frequency  coincided  with  the  required  frequency  was  indeed 
available.  The  underlined  Born  diameter  (1,030  i.m)  was  a  good  size  esti¬ 
mate,  but  only  because  the  deviation  of  the  F  estimate  from  the  effective 
transducer  center  frequency  was  acceptably  low  (2.0  vs.  2.2  MHz). 

Flaw  03  is  like  flaw  01  in  that  the  nominal  flaw  size 
(600  um)  was  properly  bracketed  by  the  underlined  estimates  (510  un  and 
990  ijm).  Flaw  04  is  also  like  flaw  01.  Flaw  05  was  underestimated  because 
a  transducer  with  an  effective  center  frequency  lower  than  2.2  MHz  was  not 
available . 


The  results  obtained  for  flaws  06  through  10  were  similar  to 
those  obtained  for  flaws  01  through  05.  Note,  however,  that  the  diameter 
of  flaw  10  was  not  properly  bracketed  by  the  underlined  estimates.  This 
casts  a  shadow  on  the  reliability  of  the  SST  to  serve  as  a  transducer 
selection  protocol. 
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Tungsten  Carbide  Inclusions  In  Titanium  Specimen 


I 


6.3.2 


The  effective  transducer  center  frequencies  required  by  the 
five  WC  inclusions  in  titanium  specimen  No.  FML  97416  are  reproduced  from 
Table  11  for  ready  reference  in  the  second  column  of  Table  14.  The  nominal 
flaw  diameters  are  listed  in  the  last  column.  Flaws  11  and  14  are  like 
flaw  01  in  that  the  nominal  flaw  sizes  were  properly  bracketed  by  the 
underlined  estimates.  Flaws  12  and  13  were  underestimated  because  a  trans¬ 
ducer  with  an  effective  center  frequency  lower  than  2.2  MHz  was  not  avail¬ 
able.  Flaw  15  is  like  flaw  10  in  that  the  flaw  diameter  was  not  properly 
bracketed  by  the  underlined  estimates. 

6.3.3  Sapphire  Inclusions  in  IN100  Specimens 


The  effective  transducer  center  frequencies  required  by  the 
eleven  sapphire  Inclusions  that  were  detected  in  IN100  specimens  Nos. 

FML  97474-1  through  -5  are  listed  in  the  second  column  of  Table  15.  The 
nominal  flaw  diameters  are  listed  in  the  last  column.  The  same  type  of 
information  is  included  in  Table  16  for  the  ten  inclusions  in  IN100  speci¬ 
men  Nos.  FML  97474-6  through  -8.  The  transducer  with  an  effective  center 
frequency  of  7  MHz  was  "discovered"  after  all  other  high-frequency  immer¬ 
sion  transducers  failed  to  detect  the  vast  majority  of  the  sapphire  inclu¬ 
sions.  All  four  flaws  that  even  this  efficient  transducer  could  not  find 
(l.e.,  flaws  17,  18,  28,  and  30)  were  small  (410  ,.m). 

Flaw  16  is  one  of  the  sixteen  sapphire  inclusions  that  only 
the  7-MHz  Immersion  LW  transducer  could  detect.  Because  the  SST  broke 
down  for  flaws  16,  26,  32,  34,  and  37  (see  Figures  A-15  through  A-17,  the 
required  frequencies  listed  in  Tables  15  and  16  for  these  410-um  diameter 
flaws  were  questionable.  The  Born  diameters  for  these  flaws  were,  there¬ 
fore,  underlined  with  caution. 


Table  14 

BORN  DIAMETER  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  TITANIUM 
SPECIMEN  NO.  FML  97416  WITH  FOUR  LONGITUDINAL-WAVE  TRANSDUCERS 


F  law 

Required 

Born  Diameter ,  _.._m 

(1) 

Nominal 

ID  No 

•  Freq.,  MHz 

2  MHz 

4  MHz 

6  MHz 

8  MHz 

Diam.  ,  .. 

1 1 

2.5 

960 

580 

490 

230 

900 

L  2 

1.8 

900 

570 

410 

250 

1,210 

1  3 

1.9 

910 

610 

390 

200 

1,210 

14 

2.5 

900 

540 

420 

2  50 

900 

15 

4.3 

o 

<r 

CN 

410 

340 

300 

480 

*-Size 

estimates  obtained 

with  the 

most  a 

ppropr iat 

e  transducers 

are 

under  1 i ned  . 


Table  15 


BORN  DIAMETER  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  IN100  SPECIMEN 
NOS.  FML  97474-1  THROUGH  -5  WITH  THREE  LONGITUDINAL-WAVE  TRANSDUCERS 


Flaw 
ID  No. 

Required 

Freq . ,  MHz 

Born 
2  MHz 

Diameter , 

4  MHz 

md) 

7  "MHz' 

Nominal 
Dlam. ,  l. 

16 

6.8 

—  (2) 

260 

410 

19 

3.2 

— 

— 

270 

790 

20 

3.2 

— 

— 

290 

790 

21 

3.1 

— 

— 

310 

790 

22 

2.2 

1,030 

450 

280 

990 

23 

2.3 

1,160 

510 

410 

990 

24 

2.2 

1,230 

460 

280 

990 

25 

3.3 

— 

430 

280 

790 

26 

6.3 

— 

— 

310 

410 

27 

2.5 

— 

480 

320 

990 

29 

3.3 

— 

— 

270 

790 

^Size  estimates  obtained 
underlined . 

2 

Ineffeclent  transducer 

with  the  most  appropriate  transducers 

(no  signal  from  flaw). 

are 

Table  16 


BORN  DIAMETER 

ESTIMATES  OBTAINED 

FOR 

THE 

INCLUSIONS  IN 

IN100  SPECIMEN 

NOS.  FML  97474 

-6  THROUGH  -8  WITH 

A  7 

MHz 

LONGITUDINAL-WAVE  TRANSDUCER 

Flaw 

Required 

Diameter 

,  ..m 

ID  No. 

Freq.,  MHz 

Born^ 1 ) 

Nominal 

31 

2.0 

390 

990 

32 

6.3 

280 

410 

33 

2.3 

300 

990 

34 

3.7 

290 

410 

35 

3.4 

280 

790 

36 

2.7 

290 

990 

37 

7.0 

260 

410 

38 

3.3 

270 

790 

39 

2.5 

300 

990 

40 

3.3 

280 

790 

^Size  estimates 

obtained  with  an  SST-a 

pproved  transducer 

are  underlined 

5  3 


Flaw  19  is  one  of  the  seven  790-.. m  diameter  inclusions  that 
only  the  7-MHz  transducer  could  detect.*  In  spite  of  this,  the  SST  did 
not  break  down  for  flaws  19,  20,  21,  29,  35,  38,  and  40  (see  Figures  A-15 
through  A-18)  so  that  the  required  frequencies  listed  in  Tables  15  and  16 
for  these  flaws  were  considered  accurate.  The  same  is,  of  course,  true 
for  flaw  25  (see  Figure  27). 

Flaw  22  is  one  of  the  three  990-..m  diameter  inclusions  that 
all  three  transducers  could  detect.  These  flaws  are  like  many  voids  and 
WC  inclusions  in  that  the  flaw  sizes  were  properly  bracketed  by  the  under¬ 
lined  Born  diameters. 

Flaw  27  is  one  of  the  five  990-.m  diameter  inclusions  that 
only  the  7-MHz  transducer  could  find.  Nevertheless,  the  SST  did  not  break 
down  for  flaws  27,  31,  33,  36,  and  39  (see  Figures  A-16  through  A-18)  so 
that  the  required  frequencies  listed  in  Tables  15  and  16  for  these  flaws 
were  considered  accurate. 

6.3.4  Real  Flaws  in  IN100  Specimens 

The  effective  transducer  center  frequencies  required  by  the 
six  real  flaws  in  INI  00  specimen  Nos.  FML  97600-1  and  -2  are  listed  in  the 
second  column  of  Table  17.  These  SST  estimates  were  suspect,  however, 
because  only  the  7-MHz  transducer  could  find  these  flaws  and  the  resulting 
magnitude  spectra  were  unmodulated  (see  Figure  A-19).  The  Born  diameters 
were,  therefore,  underlined  with  caution. 


Table  17 


BORN 

DIAMETKR  ESTIMATES 

OBTAINED  FOR 

THE  REAL 

FLAWS  IN 

IN 100  SPECIMEN 

NOS. 

FML  97600-1  and  -2 

WITH  A  7  MHz 

LONGITUDINAL-WAVE 

TRANSDUCER 

Flaw 

Requ i red 

(1) 

Actual  Flaw 

ID  No. 

Freq . ,  MHz 

Born 

Diam. ,  m 

Dimensions, 

41 

8.5 

270 

140  x  610  x  110 

42 

7.3 

250 

150  x  230  x  150 

43 

7.  3 

260 

—  (2) 

44 

7.1 

280 

370  x  500  x  90 

4  5 

8.  3 

280 

— 

46 

7.2 

260 

— 

^Size  estimates  obtained  with  the  SST-approved  transducer  are  underlined. 
^F law  has  not  been  wet  a  1 1 urg lea  1 ly  sectioned. 


*0nly  one  other  flaw  0  this  size  (flaw  25)  cou 
another  transducer. 


Id  also  be  found  with 


None  of  the  characteristics  of  these  flaws  were  known  (28) . 
For  this  reason,  three  randomly  selected  flaws  (41,  42,  and  44)  were  metal- 
lographical ly  sectioned.  No  less  than  twelve  nearly  equally  spaced  sec¬ 
tions  were  made  on  each  flaw  to  determine  its  shape,  size,  and  orientation. 
The  maximum  linear  dimensions  found  on  sequential  polishing  are  recorded 
in  the  last  column  of  Table  17.  The  underlined  Born  diameters  appear  to 
be  bracketed  by  the  minimum  and  maximum  flaw  dimensions.  Flaw  orientations 
and  dimensions  in  the  line  of  sight  had  no  effect  on  the  Born  diameters. 

The  metallographic  sectioning  procedures  were  directed  not 
only  at  verifying  the  presence  of  the  flaws  and  determining  their  shape, 
size,  and  orientation,  but  also  at  determining  the  flaw  compositions. 
Scanning  electron  microscopy  (SEM)  analysis  revealed  magnesium-rich  regions 
for  flaw  41  and  aluminum-rich  regions  for  flaw  42.  Flaw  44  was  determined 
to  be  a  void.  Flaw  compositions  had  no  apparent  effect  on  the  Born 
d iameters . 

6.4  SPOT  Size  Estimates 


The  satellite-pulse  observation  technique  was  applied  to  the  travel¬ 
ing  composite  waveforms  (i.e.,  the  doublet  trains)  obtained  by  scanning  the 
specimens  in  the  area  of  the  detected  flaw  with  high-frequency  immersion 
(or  contact)  LW  (or  SW)  transducers.  The  pulse  separation  (a)  measure¬ 
ments  were  made  directly  from  the  oscilloscope  traces  by  examiners  trained 
in  the  principles  and  application  of  the  amplitude-  and  phase-comparison 
techniques  and  the  SPOT.  The  examiners'  most  difficult  tasks  were  to 
select  the  most  appropriate  physical  model  for  the  flaw  being  examined  and 
to  discriminate  against  ghost  satellites  (material  noise  and  bondplane 
echoes).  The  selection  of  the  proper  physical  model  prior  to  the  inter¬ 
pretation  of  o  measurements  in  terms  of  flaw  dimensions  is  analogous  to 
the  selection  of  the  proper  transducer  prior  to  the  acquisition  of  the 
data  for  the  BIT.  The  SPOT  size  estimates  are  presented  separately  for 
voids,  WC  inclusions,  sapphire  inclusions,  and  real  flaws. 

6.4.1  Voids  in  Titanium  and  IN100  Specimens 


Shear-wave  amplitude  and  phase  measurements  confirmed  the 
conclusions  drawn  from  l.W  measurements  in  Paragraph  6.2.1  that  the  flaws 
in  specimen  Nos.  FML  97413  and  97414  are  spherical  voids.  A  complete  set 
of  typical  raw  waveforms  for  flaws  01  through  10  is  included  in  Appendix  C 
(see  Figures  C-l  through  C-4).  The  LW  composite  waveforms  ("signatures") 
for  flaws  03  and  08  are  reproduced  from  Figures  C-l  and  C-3,  respectively, 
for  ready  reference  in  the  upper  part  of  Figure  28.  The  SW  signatures  of 
these  flaws  are  reproduced  from  Figures  C-2  and  C-4  for  ready  reference  in 
the  upper  part  of  Figure  29.  The  heavy  trace  in  each  signature  marks  the 
beginning  and  the  end  of  the  pulse  train  (i.e.,  the  R  and  the  Sg  pulses) 
against,  the  background  of  Irrelevant  indications  (light  trace).  The  LW 
train  for  flaw  03  (dn  =  600  .  m)*  was  285  ns  long  which  translates  into  a 
size  estimate  of  620  m  [see  Equation  (7)].  The  SW  train  for  this  flaw 
was  620  ns  long,  giving  710  ..m  for  the  size  estimate  [see  Equation  (8)]. 

♦Nominal  diameter. 
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WAVEFORM 


INCLUSION  11 

c .  =  275  ns 

o  =  425  ns 
e 


-  900  ym 
d  =  910  um 
d  =  930  ym 


INCLUSION  25  -  790  ym 

r  =  360  ns  d  =  790  ym 
e 


"LAW  44  -  370  x  500  x  90  ym 

=  200  ns  d  =  440  ym 
e 


Figure  28.  Longitudinal-Wave  Signatures  of  Selected  Test  Flaws 


1 


INCLUSION  11  -  900  ym 


rLAW  44  -  370  x  500  x  90  urn 

7  =  320  ns  d  =  370  um 
e 


Figure  29.  Shear-Wave  Signatures  of  Selected  Test  Flaws 


% 


Flaw  08  (dn  =  320  .  m)  was  difficult  to  size  due  to  the  presence  of 
an  unexpected,  yet  unidentified  flaw  In  the  vicinity  of  the  Intentionally 
placed  void-  The  unexpected  asynchronous  dynamic  curves  for  the  first  and 
last  echoes  In  the  pulse  trains  (15)  gave  evidence  of  the  unintentional 
flaw  (see  the  extra  positive  half-cycle  at  the  beginning  of  the  LW  train 
in  the  second  trace  of  Figure  28).  Shear  waves  were  more  successful  than 
longitudinal  waves  in  isolating  and  subsequently  sizing  the  void  (see 
second  trace  In  Figure  29).  The  LW  and  SW  size  estimates  were  350  um  and 
280  um,  respectively. 

Typical  LW  and  SW  signatures  for  the  other  eight  voids  are 
shown  in  Figures  C-l  through  C-4.  The  averages  of  many  size  estimates 
obtained  for  the  same  void  by  varying  the  frequency,  wave  mode,  beam 
angle,  and  observation  direction  are  listed  in  Table  18. 

6.4.2  Tungsten  Carbide  Inclusions  in  Titanium  Specimen 

Shear-wave  amplitude  and  phase  measurements  confirmed  the 
flaws  in  titanium  specimen  No.  FML  97416  to  be  spherical  WC  inclusions.  A 
complete  set  of  raw  LW  and  SW  signatures  for  flaws  11  through  15  is 
included  in  Figures  C-5  and  C-6.  The  LW  and  SW  signatures  for  inclusion 
11  are  included  in  Figures  28  and  29,  respectively,  for  ready  reference. 
The  LW  train  attests  to  the  simultaneous  production  of  external  and  inter¬ 
nal  satellite  pulses  (Se  and  S^,  respectively).  The  "internal"  LW  train 


Table  18 

SIZE  ESTIMATES  OBTAINED  FOR  THE  VOIDS  IN  TITANIUM  SPECIMEN  NO. 
FML  97413  and  IN100  SPECIMEN  NO.  FML  97414  BY  THE  SATELLITE- 
PULSE  OBSERVATION  TECHNIQUE  (SPOT) 


SPOT 

Diameter,  ,.m 

Flaw 

External 

External 

Nominal 

ID  No  . 

LW  Satellites^) 

SW  Satellites^) 

Dlam.  ,  i 

01 

4  60 

370 

320 

02 

980 

1,030 

1,000 

03 

790 

700 

600 

04 

R  50 

750 

780 

05 

1,180 

1,410 

1,230 

06 

6  40 

580 

600 

07 

1,4  50 

1,470 

1,230 

08 

350 

280 

320 

09 

860 

770 

780 

10 

1,  100 

1 ,160 

1,000 

Whe  velocity  of  external  longitudinal  waves  circumventing  the  voids  was 

determined  to  be  5.3  km/sec. 
o 

The  velocity  of  external  shear  waves  circumventing  the  voids  was  deter¬ 
mined  to  be  2.9  km/sec. 


A 


n 


i 


*1 


i. 
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for  flaw  11  (dn  =  900  ..m)  was  275  ns  long,  giving  910  ..m  for  the  size  esti¬ 
mate  [see  Equation  (9)].  The  "external"  LW  train  for  this  flaw  was  425  ns 
long,  giving  930  urn  for  the  size  estimate  [see  Equation  (7)].  The  "inter¬ 
nal"  SW  train  for  this  flaw  was  625  ns  long,  giving  910/um  for  the  size 
estimate  [see  Equation  (10]. 

Typical  LW  and  SW  signatures  for  the  other  four  WC  inclu¬ 
sions  are  shown  in  Figures  C-5  and  C-6.  The  averages  of  many  size  esti¬ 
mates  obtained  for  the  same  inclusion  by  varying  the  frequency,  wave  mode, 
beam  angle,  and  observation  direction  are  listed  in  Table  19. 

6.4.3  Sapphire  Inclusions  in  IN1Q0  Specimens 

Shear-wave  amplitude  and  phase  measurements  identified  the 
flaws  in  IN100  specimen  Nos.  FML  97474-1  through  -8  to  be  spherical  sap¬ 
phire  inclusions.  A  complete  set  of  raw  LW  and  SW  signatures  for  flaws  16 
through  40  is  included  in  Figures  C-7  through  C-13.  The  LW  and  SW  signa¬ 
tures  for  inclusion  25  are  included  in  Figures  28  and  29,  respectively, 
for  ready  reference.  The  LW  train  for  flaw  25  (dn  =  790  urn)  was  360  ns 
long,  giving  790  m  for  the  size  estimate  [see  Equation  (7)].  The  SW 
train  for  this  flaw  was  685  ns  long,  giving  790  urn  for  the  size  estimate 
[see  Equation  (8)]. 

Typical  LW  and  SW  signatures  for  the  other  twenty  sapphire 
Inclusions  are  shown  in  Figures  C-7  through  C-13.  The  averages  of  many 
size  estimates  obtained  for  the  same  Inclusion  by  varying  the  frequency, 
wave  mode,  beam  angle,  and  observation  direction  are  listed  in  Tables  20 
and  21. 


Table  19 

SIZE  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  TITANIUM  SPECIMEN 

NO.  FML  97416  BY  THE  SPOT 


SPOT  Diameter,  m 


Flaw 

ID  No. 

Internal 

LW  Satellites^) 

External 

LW  Satellites^2) 

Internal 

SW  Satellites^) 

Nominal 
Diam. ,  t. 

11 

920 

930 

920 

900 

12 

1,160 

1,180 

1,150 

1,210 

13 

1,170 

1,160 

1 , 140 

1,210 

14 

910 

940 

930 

900 

1 5 

520 

480 

530 

480 

^■The  velocity  of  internal  Longitudinal  waves  in  the  tungsten  carbide 
inclusion  was  assumed  to  be  6.6  km/sec  (a  published  value). 

2The  velocity  of  external  longitudinal  waves  circumventing  the  WC  inclu¬ 
sions  was  determined  to  be  5.3  km/sec. 

The  velocity  of  internal  shear  waves  in  the  WC  inclusions  was  determined 
to  be  2.9  km/ sec . 
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Table  20 


SIZE 

ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  IN100 
NOS.  FML  97474-1  THROUGH  -5  BY  THE  SPOT 

SPOT  Diameter,  ^m 

SPECIMEN 

Flaw 

ID  No. 

External 

LW  Satellites^1) 

SW 

External 

Satellites^) 

Nominal 
Diam. ,  „m 

16 

420 

—(3) 

410 

19 

800 

— 

790 

20 

800 

— 

790 

21 

790 

— 

790 

22 

1,170 

1,210 

990 

23 

1 ,070 

1,120 

990 

24 

1,070 

1,160 

990 

25 

790 

790 

790 

26 

410 

— 

410 

27 

990 

1,080 

990 

29 

800 

760 

790 

^The  velocity  of  external  longitudinal  waves  circumvent 
Inclusions  was  assumed  to  be  5.3  km/sec. 

9 

‘'The  velocity  of  external  shear  waves  circumventing  the 


3 


slons  was  assumed  to  he  2.9  km/sec. 
No  signal  from  Inclusion. 


ing  the  sapphire 
sapphire  inclu- 


Table  21 

SIZE  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  IN100 
SPECIMEN  NOS.  FML  97474-6  THROUGH  -8  BY  THE  SPOT 


SPOT  Diameter,  ..  m 


Flaw 

Externa 1 

External 

Nominal 

ID  No- 

LW  Satellites^1) 

SW  Satellites^) 

Diam.  ,  .. 

31 

990 

1,010 

990 

32 

410 

220 

410 

33 

l  ,010 

930 

990 

34 

4  1  0 

380 

410 

35 

800 

830 

790 

36 

990 

1,0  30 

990 

37 

410 

210 

410 

38 

780 

—  (3) 

790 

39 

1  ,010 

980 

990 

40 

7  70 

7  20 

790 

1The  velocity  of  external  longitudinal  waves  circumvent: 

inclusions  was  assumed  to  be  5.3  km/sec. 

^The  velocity  of  external  shear  waves  circumventing  the 


sions  was  assumed  to  be  2.9  km/sec. 
No  signal  from  Inclusion. 


ng  the  sapphire 
sapphire  inclu- 
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6.4.4  Real  Flaws  in  IN100  Specimens 

A  complete  set  of  raw  LW  and  SW  signatures  for  flaws  41 
through  46  is  included  in  Figures  C-14  and  C-15.  The  LW  and  SW  signatures 
for  flaw  44  are  included  in  Figures  28  and  29,  respectively,  for  ready 
reference.  The  LW  train  for  flaw  44  (maximum  linear  dimensions  =  370  x 

ft 

500  x  90  .,m)  was  200  ns  long,  giving  440  .  m  for  the  equivalent  diameter 
[see  Equation  (7)].  The  SW  train  for  this  flaw  was  320  ns  long,  giving 
370  i.m  for  the  equivalent  diameter  [see  Equation  (8)]. 

Typical  LW  and  SW  signatures  for  the  other  five  real  flaws 
are  shown  in  Figures  C-14  and  C-15.  The  averages  of  many  size  estimates 
obtained  for  the  same  flaw  by  varying  the  frequency,  wave  mode,  beam 
angle,  and  observation  direction  are  listed  in  Table  22. 

6 • 5  Semlautomated  Satellite  Analysis  Procedures  (SSAP)  Size  Estimates 

The  SPOT  relies  heavily  on  the  expertise  of  the  ultrasonic  examiner 
to  extract  signal  parameters  from  an  A-scan  display.  The  satellite  pulses 
are  usually  of  low  amplitude  and  are  often  hard  to  detect  due  to  material 
noise.  The  traces  [raw  flaw  waveforms  y(t)]  in  Figure  30  reveal  the  diffi¬ 
culties  facing  the  examiner  in  recognizing  the  associated  pulses  in  y(t) 
and  Interpreting  their  separations  '  from  the  R  pulse  in  terms  of  flaw 
size  "on  the  spot."  Even  a  highly  trained  SPOT  examiner  has  difficulties 


Table  22 

SIZE  ESTIMATES  OBTAINED  FOR  THE  REAL  FLAWS  IN  IN100 
SPECIMEN  NOS.  FML  97600-1  AND  -2  BY  THE  SPOT 


SPOT  Diameter,  ..m 


Flaw 

Externa  1 

Externa  1 

Actual  Flaw 

ID  No. 

LW  Satellites^) 

SW  Satellites'1^) 

Dimensions,  nm 

41 

380 

410 

140  x 

610  x  110 

42 

390 

380 

150  x 

230  x  150 

43 

390 

320 

—  (3) 

44 

440 

390 

370  x 

500  x  90 

45 

400 

240 

— 

46 

390 

220 

— 

^The  velocity  of  external  longitudinal  waves  circumventing  the  flaws  was 
assumed  to  be  5.3  km/sec. 

2 

The  velocity  of  external  shear  waves  circumventing  the  flaws  was  assumed 
to  be  2.9  km/sec . 

Flaw  has  not  been  me ta 1 1 urgica 1 ly  sectioned. 

*The  irregular  flaws  were  modeled  as  strongly  scattering  spheres  and  sized 
in  terms  of  equivalent  diameters. 
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in  determining  where  the  pulse  train  begins  and  ends  for  some  flaws.  In 
contrast,  Figure  31  reveals  the  ease  with  which  the  spectrum  periodicities 
p  and  P  and,  subsequently,  the  flaw  sizes  can  be  obtained  from  the  flaw 
profiles.  The  automatic  detection  and  interpretation  of  the  periodicities 
clearly  present  in  11(f)  is  a  significant  step  toward  eliminating  the  need 
for  a  highly  trained  examiner  in  the  inspection  loop  ( 29 ).  The  semiauto- 
mated  satellite  analysis  procedures  (SSAP)  detailed  in  Table  23  are  the 
first  step  in  automating  the  SPOT  to  make  its  use  practical  by  less  skilled 
ultrasonic  examiners.  Better  pulse  resolution  can  be  achieved  not  only  by 
using  shear  as  well  as  longitudinal  waves,  but  also  by  deconvolving  y(t) 
with  r(t)  (i.e.,  its  own  front-surface  echo).  Not  only  can  the  spectrum 
periodicities  be  read  easier  than  the  pulse  separations,  but  the  accuracy 
of  p/P  measurements  is  also  improved  over  that  of  a  measurements. 

In  generating  the  data  shown  in  Tables  24  through  26  for  the  32 
detectable  inclusions,  the  SW  measurements  were  used  to  confirm  the 
physical  model  applied  to  and  the  size  estimates  derived  from  the  LW  mea¬ 
surements.  A  complete  set  of  profiles  for  these  inclusions  is  contained 
in  Appendix  D. 


Table  23 

EXPERIMENTAL  PROTOCOL  FOR  THE  SEMI AUTOMATED  SATELLITE  ANALYSIS 

PROCEDURES  ( SSAP) 


_ Description _ 

Record  flaw  waveform  y(t)  with  any  high-frequency  LW 
transducer 

Digitize  and  store  y(t) 

Gate  out  the  R  pulse  of  y(t)  to  serve  as  reference 
waveform  r(t) 

Digitize  and  store  r(t) 

Deconvolve  y(t)  with  r(t) 

Determine  periodicity  p  of  magnitude  spectrum  11(f)  in  MHz 

Classify  flaw  as  strong  (void  or  sapphire)  or  weak  (WC) 

Estimate  flaw  size  in  .  m  by  2,'80/p  (strong  flaw)  or 
3,300/p  (weak  flaw) 

Repeat  steps  1  through  7  with  any  high-frequency  SW 
t  ransducer 


Estimate  flaw  size  in  m  by  1,130/1’  (strong  flaw)  or 
1,4  30/1’  (weak  flaw) 


L  WAVE  PROBE 


S  WAVE  PROBE 


2 


3 


7 


7 


1  2 


ligure  il.  1,011)’  i  t  ml  i  na  1  -  and  Shear-Wave  Profiles  of  Tungsten 
Carbide  Inclusions  11,  12,  and  IS 


Table  24 


SIZE  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  TITANIUM  SPECIMEN  NO. 

FML  97416  BY  THE  SSAP 


Flaw 

SSAP 

Diameter,  ;.m 

Nominal 

ID  No. 

L  Waves U ) 

S  Waves(^) 

Diam . ,  „  m 

11 

940 

910 

900 

12 

L  ,270 

1,120 

1,210 

13 

1,220 

1,120 

1,210 

14 

940 

910 

900 

15 

450 

520 

480 

^High-frequc  icy  contact  transducer. 
^High-f requency  immersion  transducer. 


Table  25 

SIZE  ESTIMATES  OBTAINED  FOR  THE  INCLUSIONS  IN  IN100  SPECIMEN 
NOS.  FML  97474-1  THROUGH  -8  BY  THE  SSAP 


F 1  aw 

Diameter,  .m 

Flaw 

Diameter 

,  ..m 

ID  No. 

SSAP^ 1 ) 

Nominal 

ID  No. 

SSAPO) 

Nominal 

16 

520 

410 

31 

1,100 

990 

19 

690 

790 

32 

470 

410 

20 

o70 

790 

33 

1  ,  100 

990 

21 

650 

790 

34 

510 

410 

22 

1  , 1 00 

990 

35 

730 

790 

23 

1 , 0  50 

990 

3  b 

1,050 

990 

24 

960 

990 

37 

440 

410 

25 

690 

790 

38 

690 

790 

2b 

500 

410 

39 

1,100 

990 

27 

920 

990 

40 

710 

790 

29 

710 

790 

^High-frequency  contact  LW  transducer. 


Table  26 


SIZE  ESTIMATES  OBTAINED  FOR  THE  REAL  FLAWS  IN  IN100  SPECIMEN 
NOS.  FML  97600-1  AND  -2  BY  THE  SSAP 


Flaw 

ID  No. 

SSAP  Diam. ,  .  m^1) 

Actual  Flaw 
Dimensions , 

41 

230  -  450 

140  x  610  x 

110 

42 

220  -  430 

150  x  230  x 

150 

43 

250  -  500 

—  (2) 

44 

260  -  520 

370  x  500  x 

90 

45 

250  -  500 

— 

46 

240  -  480 

— 

^High-frequency  contact  LW  transducer. 

^Flaw  has  not  been  me ta 1 lurglcal ly  sectioned. 
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7.  DISCUSSION 


7 . L  General 

In  this  section,  the  size  estimates  obtained  under  blind  test  con¬ 
ditions  for  the  42  detected  flaws  by  the  primary  and  the  additional  flaw 
characterization  procedures  are  discussed.  The  Born  inversion  procedures 
and  those  based  on  the  detection  of  associated  pulses  (i.e.,  the  SPOT  and 
SSAP)  differ  fundamentally  because  the  former  procedures  attempt  to  deter¬ 
mine  the  flaw  size  solely  from  the  parameters  of  the  front-surface  echo 
(1*22,24). 

7 . 2  Born  Inversion  Procedures 

As  shown  in  Table  4,  the  BIT  needs  supplementation  by  the  SST  which, 
in  turn,  uses  the  results  of  the  amplitude-  and  phase-comparison  techniques 
as  well  as  those  of  deconvolution  in  selecting  the  most  appropriate  trans¬ 
ducers)  to  use  in  the  acquisition  of  waveforms  for  Born  inversion.  The 
amplitude-  and  phase-comparison  techniques  were  successfully  applied  to 
the  raw  flaw  waveforms  obtained  with  the  high-frequency  LW  transducers. 

All  36  detectable  flaws  in  the  diffusion-bonded  test  specimens  were  cor¬ 
rectly  identified  as  spherical  voids,  WC  inclusions,  or  sapphire  inclusions. 
The  appropriateness  of  an  LW  transducer  for  inclusion  into  the  Born  mea¬ 
surement  protocol  was  determined  by  preliminary  examination  of  the  resulting 
deconvolved  flaw  spectrum  (magnitude  spectrum  11(f)] .  A  transducer  whose 
spectral  response  W(f)  encompassed  more  than  one  of  the  flaw's  resonance 
peaks  was  declared  inappropriate  for  acquiring  Born  inversion  data. 

Flaw  02  serves  to  illustrate  this  transducer  selection  protocol.  All  of 
the  high-frequency  probes  used  to  interrogate  this  flaw  were  considered 
inappropriate  (see  Figure  22).  Note  that  the  bandwidth  of  the  selected 
2-MHz  transducer  is  well  centered  on  the  first  peak  of  H(f). 

When  none  of  the  available  transducers'  bandwidth  was  centered  on 
the  flaw's  fundamental  (n=l)  resonance  frequency,  two  transducers  were  con¬ 
sidered  appropriate  for  acquiring  Born  inversion  data.  Flaws  01  and  03 
serve  to  illustrate  this  transducer  selection  protocol.  Both  the  4-MHz 
and  8 -MHz  transducers*  were  considered  appropriate  for  flaw  01  because  the 
Born  frequency  was  estimated  to  be  6  MHz.  A  Born  frequency  estimate  of 
2.8  MHz  influenced  the  selection  of  the  2-MHz  and  4-MHz  transducers  for 
flaw  03.  Because  a  large  error  in  the  F  estimate  by  the  SST  could  have 
led  to  an  unacceptable  error  in  the  D  estimate  by  the  BIT,  Born  Inversions 
were  performed  for  the  waveforms  acquired  for  a  given  flaw  by  all  three 
standard-frequency  transducers.  However,  only  those  Born  diameters  that 
were  obtained  with  the  SST-approved  transducers  are  underlined  in 
Tables  25  hrough  29. 

The  BIT  yielded  size  estimates  (D)  with  large  inaccuracy  for  the  15 
voids  and  WC  inclusions  as  well  as  nine  of  the  sapphire  inclusions  for 

*Tho  effective  center  frequencies  of  the  available  2.25,  5,  and  10  MHz 
transducers  were  approx  I  mate ly  2,  4,  and  8  MHz,  respectively. 
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which  at  least  one  transducer  with  an  effective  center  frequency  near  the 
flaw's  fundamental  resonance  was  available  (coincidence  rule).  Figure  32 
shows  that  the  correlation  of  the  Born  diameters  with  the  nominal  flaw 
sizes  is  barely  significant  (correlation  coefficient  0.81).  The  sizing 
errors  for  the  24  flaws  were  such  that  310  m  (L2  mils)  need  to  be  added 
to  D  in  order  to  be  almost  certain  (i.e.,  95  percent  sure)  that  the  actual 
flaw  size  does  not  exceed  this  conservative  size  estimate. 

The  Born  diameters  were  determined  to  be  independent  of  flaw  com¬ 
position,  size,  and  shape,  but  linearly  related  to  the  period  of  the 
ultrasonic  wave  (T)  according  to  the  equation 

D  =  1.95  T  (11) 

where  D  is  in  micrometers  and  T  is  in  nanoseconds.  This  relationship  was 
found  empirically  by  plotting  the  averages  of  the  37  Born  diameters  under¬ 
lined  in  Tables  13  through  17  against  the  effective  wave  periods  (see 
Figure  33).  Equation  (LI)  can  be  better  understood  following  the  deriva¬ 
tion  of  the  coincidence  rule  formulated  in  the  time  domain.  It  can  be 
shown  (see  Table  27)  that  the  bandwidth  requirement  for  the  Born  trans¬ 
ducer  (i.e.,  0.5<ka<2.5)  reduces  to 

T  =  og  (12) 

Matching  the  transducer  center  frequency  to  the  flaw's  fundamental  reso¬ 
nance  frequency  thus  guarantees  the  smooth  blending  of  the  external 
satellite  pulse,  if  one  is  produced,  into  the  natural  ringdown  (decay)  of 
the  front-surface  echo  (R  pulse).  This  explains  why  the  presence  of  the 
satellite  pulse  of  a  strongly  scattering  flaw  is  not  apparent  when  the 
flaw  waveform  is  acquired  with  the  "proper"  transducer.  Since  the  sepa¬ 
ration  of  the  third  half-cycle  from  the  first  half-cycle  (R  pulse)  in  the 

flaw  waveform  acquired  by  the  Born  transducer  is  equal  to  T,  the  BIA 
treats  the  third  half-cycle  in  y(t)  as  if  it  were  the  external  satellite 

pulse.  The  algorithm  does  not  actually  need  a  satellite  pulse  as  long  as 

the  ringdown  of  the  R  pulse  lasts  for  at  least  one  period  (i.e.,  for  two 
more  half-cycles). 


Table  27 

DERIVATION  OF  THE  TIME-DOMAIN  COINCIDENCE  RULE  FOR  A  STRONG  SCATTERER 


Quantity  _ Equation _ _ Remarks 


Born  bandwidth 

0 .  'xka-;  2 . 5 

Transducer  selection 

Born  wave  number 

Ka  =  (0. 5x2.5)l/2  =1.1 

Geometric  mean 

Born  frequency 

F  =  1 .  1  c  /  *  d 

f  =  kc/2- 

Born  period 

T  =  2.8  d/c 

T  =  L/F 

Pulse  separation 

•  ’  2.8  d/c 

Equation  (6) 

Coincidence  rule 

T  =  ’ 

c 

Creeping  wave  reinforce: 
fleeted  wave  (construct 
interference) 

Figure  32.  Correlation  of  Nominal  Flaw  Sizes  With  Estimates 
Obtained  by  the  Born  Inversion  Procedures  (BIP) 


Born  diameters  were  obtained  with  acceptable  inaccuracy  only  when  all  of 
these  conditions  were  met-  The  main  difficulty  with  the  low-frequency 
(ka-0.5)  backscattering  measurements  is  that  flaws  in  materials  exhibiting 
significant  grain  scattering  can  be  easily  missed.  Four  of  the  smallest 
sapphire  inclusions  in  the  IN100  test  specimens  were  not  found  at  all. 

The  BIP  completely  broke  down  for  22  of  the  1N100  flaws  because  these 
flaws  were  not  detectable  by  more  than  one  transducer.  The  coincidence 
rule  was  violated  for  most  of  these  flaws. 

The  Rockwell  Born  diameter  estimates  are  generally  higher  than  the 
Ames  estimates  (see  Table  6).  This  is  attributed  mainly  to  the  difference 
in  the  adapted  ZOTA.  The  "proper"  frequency  for  the  standard  flaw  (800-;. m 
diameter  void)  is  calculated  by  Equation  (1)  to  be  2.8  MHz  (see  Figure  17 
for  the  Ames  BIP  results  obtained  with  a  2.5-M!iz  transducer).  Neverthe¬ 
less,  the  Rockwell  BIP  produced  a  size  estimate  of  760  m  for  this  flaw 
with  a  4.3-MHz  "optimum"  transducer  by  resorting  to  an  empirical  ZOTA  (see 
Figure  19).  The  r  parameter  was  adjusted  until  the  flatness  of  the  charac¬ 
teristic  function  J(r)  was  maximum  at  r=0.  The  effect  was  an  increase  in 
the  Born  diameter  from  dQ  =  510  „m  to  D  =  760  i.m.*  It  was  this  "success" 
of  the  BIP  that  prompted  us  to  consider  testing  yet  another  ZOTA  (30). 

This  ZOTA  is  based  on  the  Born  model  and  it  simply  requires  the  time  shift 
t  to  be  oe/2.  The  Born  diameters  obtained  for  five  representative  inclu¬ 
sions  using  the  "correct"  (Born)  r  are  compared  with  those  obtained  with 
the  Ames  BIP  in  Table  29.  It  is  concluded  that  neither  ZOTA  is  capable  of 
assuring  the  reproducibility  of  Born  diameters  obtained  for  the  same  flaw 
with  different  transducers.  The  influence  of  the  transducer  center  fre¬ 
quency  on  the  Born  diameters,  exhibited  in  Figure  33,  is  much  stronger  and 
more  predictable  than  that  of  x  . 


Table  29 

DEPENDENCE  OF  BORN  DIAMETER  ESTIMATES  ON  METHOD 
FOR  LOCATING  THE  FLAW  CENTER 


Flaw 

Flaw 

Born 

Diameter,  .m 

Nomina  1 

ID  No. 

Type 

Ames  t 

Born  ■  (=  i e/2 ) 

Diam.  ,  i.  m 

11 

WC 

310 

310 

900 

13 

wc 

380 

730 

1,210 

15 

WC 

300 

290 

480 

36 

Sapph 1  re 

340 

700 

990 

37 

Sapph 1  re 

250 

700 

410 

*Note,  however,  that  when  the  transducer  was  "proper,"  the  same  empirical 
procedures  for  locating  the  flaw  center  increased  the  Born  diameter  from 
dQ  =  800  ;  m  (the  right  answer)  to  D  =  1,230  .  m  (see  Figure  21). 


Satellite-Pulse  Observation  Technique 


The  time-domain  formalism  is  rich  In  Insights  to  the  flaw-ultrasound 
interaction  and  permits  simple  solutions  to  the  inverse  scattering  prob¬ 
lem  (4).  The  flaw's  response  to  short  bursts  of  high-frequency  insonifica- 
tion  can  be  examined  in  conjunction  with  a  physical  model  directly  on  the 
instrument  screen  (SPOT)  or  analyzed  by  a  semiautomated  pattern  recogni¬ 
tion  algorithm  (see  Subsection  7.4  for  SSAP). 

Three  physical  models  were  developed  to  aid  Lhe  ultrasonic  examiner 
in  the  recognition  and  real-time  interpretation  of  the  satellite  pulses 
(both  internal  and  external)  obtained  for  the  voids,  WC  inclusions,  and 
sapphire  inclusions  in  the  test  specimens  with  h igh-f requency  LW  and  SW 
transducers.  The  transducer  center  frequency  per  se  is  immaterial  so  long 
as  its  transient  response  is  short  enough  to  permit  resolution  of  the  ESP. 

The  SPOT  was  supplemented  by  the  phase-comparison  technique  in  the 
flaw  identification  task.  The  polarities  of  the  associated  R  and  S  pulses 
for  the  three  flaw  types  were  as  expected,  i.e.,  the  reflected  pulse  was 
either  positive  or  negative,  the  internal  satellite  was  always  positive, 
and  the  ESP  was  always  negative.*  These  polarity  rules  were  of  immense 
value  in  both  flaw  identification  and  sizing. 

All  flaws  were  "strong"  enough  to  support  externally  scattered  waves. 
This  result  was  not  expected.  External  satellite  pulses  were  observed  in 
all  but  one  experimental  situation  (shear  waves  incident  on  WC  inclusions). 
Since  the  separation  of  the  ESP  from  the  R  pulse  is  independent  from  flaw 
composition,  the  inverse  scattering  formulas  Involving  [i.e..  Equations 
(7)  and  (8)]  are  applicable  to  all  types  of  flaws  with  smooth  boundaries. 

Physical  models  were  also  developed  for  the  production  of  bondplane 
echoes  from  the  IN100  specimens  so  that  these  irrelevant  indications  could 
be  promptly  discounted  in  search  of  the  often  weak  satellites  of  the  R 
pulse  (satellite-pulse  hunt). 

The  SPOT  diameters  listed  in  Tables  18  and  19  for  the  ten  voids  and 
five  WC  inclusions  represent  the  averages  of  many  measurements  (over  50  in 
most  cases)  performed  with  a  multitude  of  transducers  (shear  vs.  longitudi¬ 
nal,  medium  vs.  high  frequency,  Immersion  vs.  contact,  and  low  vs.  high 
beam).  In  addition,  each  flaw  was  viewed  with  the  same  transducer  from 
many  directions.  The  1,W  and  SW  ~  measurements  were  treated  on  an  equal 
footing.  The  shear  waves  were  more  effective  in  producing  Se  pulses  from 
voids,  and  the  longitudinal  waves  produced  stronger  pulses  from  WC 
inclusions.  This  mul t i probe-and-mul t Ipulse  approach  has  proven  to  be  a 
good  substitute  for  deconvolution,  since  we  capitalized  on  the  independence 
of  '  measurements  from  transducer  characteristics  (including  probe 


*There  were  four  exceptions  out  of  the  42  cases  to  the  last  polarity  rule. 


position).  The  measurements  were  repeated  in  different  sequence  to 
assure  reproducibility  of  the  results. 

The  procedure  qualification  tests  included  one  situation  in  which 
the  intended  flaw  (void  08)  was  in  the  "shadow  zone"  of  an  unintentional, 
yet  unidentified  reflector.  The  SW  measurements  were  most  helpful  in 
isolating  and  subsequently  sizing  this  flaw. 

The  SPOT  size  estimates  shown  in  Figure  34  for  the  10  voids  and  the 
26  detectable  inclusions  In  the  diffusion-bonded  test  specimens  are  seen 
to  be  in  good  agreement  with  the  flaw  dimensions  as  reported  by  the  speci¬ 
men  manufacturer.  The  correlation  of  the  SPOT  estimates  with  the  nominal 
flaw  sizes  is  significant  (correlation  coefficient  0.97).  For  this  reason, 
the  metallographic  sectioning  of  any  of  these  flaws  was  waived.  Future 
flaw  characterization  module  development  programs  will  benefit  from  the 
availability  of  these  (now  standard)  specimens.  The  sizing  errors  for  the 
36  test  flaws  were  such  that  merely  100  ..m  (4  mils)  need  to  be  added  to 
the  SPOT  estimate  in  order  to  be  almost  certain  that  the  actual  flaw  size 
does  not  exceed  this  conservative  estimate. 

Three  randomly  selected  real  flaws  (41,  42,  and  44)  were  metallo- 
graphically  sectioned  to  determine  how  well  the  SPOT  and  the  SSAP  per¬ 
formed  in  a  "real"  test.*  Real  flaws  are  complex  in  shape  and  disagree¬ 
ments  between  actual  and  estimated  dimensions  were  expected.  The  scanning 
electron  microscope  (SEM)  photographs  of  the  sectioned  flaws  are  shown  in 
Figure  35.  Hand  grinding  was  done  in  the  YZ  plane.  In  each  case,  the 
largest  flaw  dimension  was  in  the  X  direction.  The  aluminum-rich  flaw  42 
was  the  closest  to  being  spherical.  The  magnesium-rich  flaw  41  was 
needle-like.  Flaw  44  (a  void)  was  almost  two-dimensional  (planar). 

Figures  C-14  and  C- L 5  contain  the  LW  and  SW  signatures,  respec¬ 
tively,  for  the  sectioned  flaws.  The  through-wall  (depth)  and  length 
dimensions  cannot  b<  ascertained  from  the  backscattered  waveforms.  It  was 
expected,  however,  that  the  interpretation  of  extra  pulses  in  y(t)  will  be 
similar  to  that  described  in  Subsection  2.7  and  that  in  many  cases  the 
real  flaws  can  be  approximated  by  regularly  shaped  flaws. 

The  SPOT  and  SSAP  size  estimates  listed  in  Table  30  appear  to  be 
conservative  lor  flaw  42  and  bracketed  by  the  largest  two  dimensions  of 
flaws  41  and  44.  All  real  flaws  were  issuraed  to  be  strong  spherical 
scatterers.  These  results  are  good  considering  that  these  flaws  were 
neither  spherical  nor  homogeneous.  It  remains  to  be  seen  whether  the 
other  three  real  t  1  iws  are  also  Irregular  and  inhomogeneous. 

7.4  Semiaut  omat  ed  Sate  1 1 1 te  Analysis  Procedures 

It  was  not  unt  il  well  into  the  program,  when  the  size  estimates 
obtained  by  the  primary  and  additional  flaw  characterization  procedures 
for  the  first  ten  test  flaws  were  compared  to  their  nominal  dimensions, 

*Flaws  43,  45,  and  4h  were  held  tor  further  measurements. 
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FLAW  41 

140  x  610  x  100  ..m 


FLAW  42 


150  x  230  x  150  i.m 


FLAW  44 


370  x  500  x  90  ,m 


1  j  ■  rosropy  ( S KM )  Photographs 


Table  30 


COMPARISON  OK  ACTUAL  AND  ESTIMATED  SPOT  AND  SSAP 
DIMENSIONS  OBTAINED  FOR  THREE  REAL  FLAWS 


Flaw  Equivalent  Diameter,  ..m  Largest  Two  Flaw 


ID  No. 

SPOT 

SSAP 

Dimensions,  .  m 

Compos i t ion 

41 

380,  410 

230,  450 

140,  610 

Magnesium 

42 

380,  390 

220,  430 

150,  230 

Aluminum 

44 

390,  440 

260,  520 

370,  500 

Void 

that  attempts  were  made  to  automate  or  systematize  the  qualified  SPOT  pro¬ 
cedures  to  the  extent  possible  In  order  to  assure  reproducibility  and  reli¬ 
ability.  The  remaining  32  test  flaws  were,  therefore,  used  not  only  to 
further  test  the  applicability  of  the  SST,  BIT,  and  SPOT  to  inclusions  but 
also  to  quantify  the  performance  degradation  associated  with  the  elimina¬ 
tion  of  the  SST  and/or  SPOT  examiner  from  the  inspection  loop. 

The  SSAP  size  estimates  shown  in  Figure  36  for  the  five  WC  inclusions 
and  the  21  detectable  sapphire  inclusions  are  seen  to  be  less  accurate — by 
AO  percent  to  be  exact — than  those  obtained  directly  from  the  instrument 
screen  by  a  trained  examiner  (see  Figure  34).  The  correlation  of  the  SSAP 
estimates  with  the  nominal  flaw  sizes  remains  significant  (correlation 
coefficient  0.94).  The  sizing  errors  for  the  26  inclusions  were  such  that 
merely  140  .m  (5.5  mils)  need  to  be  added  to  the  SSAP  estimate  in  order  to 
be  almost  sure  that  the  actual  flaw  size  does  not  exceed  this  conservative 
estimate.  Sizing  errors  of  this  magnitude  are  acceptable  for  an  automated 
flaw  characterization  module.  The  performance  of  the  SSAP  was  degraded 
relative  to  that  of  the  SPOT  mainly  because  the  SSAP  cannot  differentiate 
against  the  spurious  echoes  from  the  bondplanes  in  the  IN100  specimens  and 
only  longitudinal  waves  were  used  to  interrogate  the  sapphire  inclusions. 

Of  course,  bond  plane  indications  were  not  a  problem  for  the  SSAP  in 
the  analysts  of  the  profiles  obtained  for  the  six  real  flaws  (see  Fig¬ 
ure  D —  7 ) .  The  SSAP  size  estimates  obtained  for  these  flaws  are  included 
In  Table  30.  As  expected,  deconvolution  with  the  front-surface  echo  pro¬ 
duced  estimates  that  wore  more  precise  (smaller)  than  the  SPOT  estimates. 
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8.  CONCLUSIONS 


The  main  program  objective  was  to  evaluate  the  supplemented  BIT  (Born 
inversion  technique)  in  the  forms  developed  at  the  Ames  Laboratory  and  the 
Rockwell  Science  Center  against  40  or  more  internal  bulk  flaws  in  titanium 
and  IN100  alloys  under  simulated  field  (blind-test)  conditions.  Due  to 
its  sole  reliance  on  the  front-surface  echo  (R  pulse),  the  BIP  (Born  inver¬ 
sion  procedures)  yield  essentially  the  same  result  for  voids  and  inclusions 
when  the  same  transducer  is  used.  The  Born  diameter  estimates  (D)  were 
determined  to  be  independent  of  not  only  flaw  composition,  but  also  shape, 
size,  and  orientation.  Instead,  the  size  estimates  were  determined  to  be 
directly  related  to  the  period  of  the  incident  ultrasonic  wave  (T).  Each 
transducer  yields  a  different  Born  diameter  according  to  the  relationship 
D  =  1.95T,  where  D  is  in  micrometers  and  T  is  in  nanoseconds. 

The  BIA  (Born  inversion  algorithm)  expects  two  pulses  of  comparable  ampli¬ 
tude  from  the  flaw  which  is  assumed  to  be  a  weak  scatterer  (Born  flaw) . 

The  algorithm  correctly  treats  the  first  half-cycle  of  the  impulse  response 
of  a  real  flaw  as  the  front-surface  echo,  but  mistakes  its  third  half-cycle 
for  a  back-surface  echo  (S^  in  Figure  C-5)  or  a  creeping-wave  return  (Se 
in  Figure  C-l).  Thus,  short-period  reflected  pulses  give  rise  to  small 
Born  diameters  and  long-period  reflected  pulses  yield  large  size  estimates. 
For  this  reason,  the  BIP  results  could  not  be  reproduced  with  different 
transducers.  Unless  there  was  a  transducer  available  for  transmitting 
ultrasonic  waves  with  a  period  equal  to  the  separation  oe  between  the 
associated  reflected  and  externally  scattered  pulses  (coincidence  rule, 

T  =  oe),  the  Born  diameter  was  unreliable.  This  means  that  of  the  trans¬ 
ducers  used  to  obtain  the  ten  flaw  waveforms  shown  in  Figures  C-l  and  C-5, 
only  the  transducer  used  on  flaw  01  was  "proper”  for  acquiring  the  Born 
inversion  data.  For  the  other  nine  flaws,  the  effective  transducer  center 
frequency  had  to  be  lowered  by  at  least  50  percent  in  order  to  obtain 
reliable  Born  diameters. 

Presorting  of  the  three  standard-frequency  transducers  for  acquiring  the 
Born  inversion  data  was  done  by  applying  the  qualitative  SST  (spectrum 
splicing  technique)  to  the  deconvolved  flaw  spectra  obtained  with  the 
available  transducers.  A  transducer  whose  bandwidth  encompassed  more  than 
one  of  the  flaw's  resonance  peaks  was  ruled  inappropriate  for  acquiring 
Born  inversion  data.  Efficient  appropriate  transducers  were  available  for 
24  of  the  36  flaws  that  were  detected  in  the  diffusion-bonded  test  speci¬ 
mens.*  The  correlation  of  the  resulting  Born  diameters  with  the  nominal 
flaw  dimensions  was  barely  significant  (correlation  coefficient  0.81). 

The  size  estimates  were  accurate  only  to  within  370  um  (see  Figure  32). 

One  certain  way  to  improve  the  reliability  and  accuracy  of  the  Born  size 
estimates  is  to  include  non-standard  as  well  as  standard  frequency  trans¬ 
ducers  in  the  measurement  protocol.  This  is  not  practical. 


*The  BIP  broke  down  completely  for  the  12  sapphire  inclusions  that  could 
not  be  detected  with  the  required  low-frequency  transducers. 
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The  sensitivities/insensitivities  of  and  the  frequencies  preferred  by  the 
five  techniques  incorporated  into  the  primary  and  additional  flaw  charac¬ 
terization  procedures  are  summarized  in  Table  31.  The  amplitude-  and 
phase-comparison  techniques  were  used  to  classify  the  flaw  as  either  strong 
or  weak.  The  SST  used  this  information  in  determining  the  most  appropriate 
transducer (s)  for  acquiring  the  Born  inversion  data  for  a  given  flaw.  The 
BIP  and  the  SPOT  (satellite-pulse  observation  technique)  use  different 
information  to  obtain  the  size  estimate.  There  is  a  great  deal  of  flaw 
information  thrown  away  by  the  BIP.  By  the  time  the  effective  transducer 
center  frequency  is  lowered  so  that  T  =  oe»  the  ESP  (external  satellite 
pulse)  can  no  longer  be  distinguished  from  the  ringdown  of  the  R  pulse. 

In  contrast  to  the  BIP,  the  SPOT  is  sensitive  to  Z  (acoustical  impedance 
of  flaw)  and  d,  but  insensitive  to  T  as  long  as  the  frequency  of  the  inci¬ 
dent  wave  is  high  (fr  >  4  MHz).  It  is  concluded  that  multiple  measure¬ 
ments  made  with  waves  of  different  mode,  period,  and  incidence  angle  and 
no  less  than  three  physical  models  for  ultrasound-flaw  interactions  (one 
for  voids  and  the  other  two  for  weakly  and  strongly  scattering  inclusions) 
are  necessary  to  solve  the  inverse  problem  for  flaw  type  and  size  in  an 
unambiguous  manner. 


Table  31 

COMPARISON  OF  THE  FIVE  ULTRASONIC  TECHNIQUES  APPLIED 
TO  THE  TEST  FLAWS 


Technique 


Sensitive 

To 


Insensitive 

To 


Preferred 
Frequency  Range 


Amplitude  comparison 
(qualitative) 

Phase  comparison 
(qualitative) 

Spectrum  splicing 
(qualitative) 

Born  inversion 

(quantitative) 

Satellite  pulse 
(quantitative) 


dU),  Z<2) 


T(3) 

d,  T 


d,  Z 


Z,  d 


Z,  d 


N.A.(4> 

N.A. 

Low  to  high 
Born  frequency 
Medium  to  high 


1-Flaw  diameter 

^Acoustical  impedance  of  flaw 
^Wave  period 
^Not  applicable 
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Only  high-fidelity  transducers  that  emphasized  the  important  physical  vari¬ 
ables  involved  in  the  flaw-ultrasound  interaction  process  were  included  in 
the  SPOT  measurement  protocol.  Space-time-frequency  analyses  of  the  infor¬ 
mation  obtained  "on  the  spot"  about  the  flaw  over  the  4  to  8  MHz  frequency 
range  resulted  in  a  significantly  better  characterization  of  the  test  flaws 
than  has  been  possible  with  the  B1P.  All  of  the  42  detected  flaws  were 
sized  by  noting  the  separation  of  the  ESP  from  the  specularly  reflected 
pulse.  Internal  satellite  pulses  in  the  case  of  the  five  WC  inclusions 
provided  confirmatory  size  estimates.  The  six  real  flaws  were  too  small 
for  determination  of  their  shape  and  orientation.  The  size  estimates  for 
these  flaws  were  conservative.  The  sizing  errors  for  the  36  flaws  in  the 
diffusion-bonded  specimens  were  such  that  merely  100  wm  need  to  be  added 
to  the  SPOT  estimate  in  order  to  be  almost  sure  that  the  actual  flaw  is 
not  larger  than  this  conservative  size  estimate. 

The  satellite  pulses  can  be  best  enhanced  by  the  deconvolution  process  if 
the  front-surface  echo  [i.e.,  the  first  three  half-cycles  in  y(t)]  is  used 
as  the  reference  waveform.  Unlike  the  results  of  doublet-separation  (o) 
measurements,  the  results  of  resonance-peak-separation  (p  or  P)  meassure- 
ments  can  be  averaged  during  specimen  scanning.  Thus,  satellite  pulses 
can  best  be  detected  by  analyzing  the  spatially  averaged  magnitude  spectra 
H(f )  for  the  degree  and  extent  of  modulation.  As  long  as  the  R-S^  doublet 
of  inclusion  11  (drawn  with  a  heavy  line  in  the  uppermost  trace  of  Fig¬ 
ure  C-6)  does  not  leave  the  gate,  its  magnitude  spectrum  (see  upper  right 
trace  in  Figure  D-l)  is  time  invariant.  The  essence  of  SSAP  (semiautomated 
satellite  analysis  procedures)  is  that  the  peak  separation  P^  can  be  deter¬ 
mined  automatically  while  the  accurate  assessment  of  the  doublet  separation 
oj  (**  1/Pj)  requires  a  judgment  by  the  examiner.  It  can  be  concluded  that 
the  semiautomated  satellite  analysis  procedures: 

Must  be  supplemented  by  techniques  to  classify  the  flaw  as  strong  or 
weak. 

Do  not  require  the  examiner  to  interpret  the  flaw  waveform  directly 
from  the  instrument  screen. 

■  Yield  more  precise  size  estimates  than  the  SPOT. 

•  Make  more  complete  use  of  the  available  flaw  information  than  the 
Born  inversion  procedures. 

Are  better  candidates  for  automated  application  than  the  Born 
inversion  procedures. 
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Figure  A-2.  Magnitude  Spectra  Obtained  for  Titanium  Void  No .  02 
with  Four  Longitudinal-Wave  Transducers 


Figure  A-4.  Magnitude  Spectra  Obtained  for  Titanium  Void  No.  04 
with  Four  Longitudinal-Wave  Transducers 
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Figure  A-19.  Magnitude  Spectra  Obtained  for  IN100  Flaws  Nos.  41 

through  46  with  a  7  MHz  Longitudinal-Wave  Transducer 
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Figure  B-l.  Input-Output  Relationships  Yielded  by  the  Born  Inversion 
Algorithm  (BIA,  Flaws  01  and  02).  Note  that  k*a*=0.92, 
0.80,  and  0.96. 
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Figure  B-2.  Input-Output  Relationships  Yielded  by  the  BIA  (Flaws  02 
and  03).  Note  that  k*a*“0.92,  0.96,  and  0.89. 
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Figure  B-4.  Input-Output  Relationships  Yielded  by  the  BIA  (Flaws  10 
and  11).  Note  that  k*a*«0.94,  0.95,  and  0.94. 
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Figure  B-6.  Input-Output  Relationships  Yielded  by  the  BIA  (Flaws  16 
and  22).  Note  that  k*a*=0.88,  0.85,  and  0.92. 
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Figure  B-7.  Input-Output  Relationships  Yielded  by  the  BIA  (Flaws  41 
through  43).  Note  that  k*a*=0.92,  0.86,  and  0.90. 
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Figure  B-8.  Input-Output  Relationships  Yielded  by  the  BIA  (Flaws 
through  46).  Note  that  k*a*=0.90,  0.88,  and  0.90. 


APPENDIX  C 

LONGITUDINAL-  AND  SHEAR-WAVE  SIGNATURES 
USED  IN  THE  APPLICATION  OF  THE 
SATELLITE-PULSE  OBSERVATION  TECHNIQUE 
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Figure  07.  Longitudinal-Wave  Signatures  of  the  Inclusions 
in  IN100  Specimen  Nos.  FML  97474-1  through  -1 
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Figure  09.  Longitudinal-Wave  Signatures  of  the  Inclusions 
in  IN100  Specimen  Nos.  FML  97474-5  through  -7 
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Figure  C-15.  Shear-Wave  Signatures  of  the  Real  Flaws  in 
IN100  Specimen  Nos.  FML  97600-1  and  -2 
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APPENDIX  D 

LONGITUDINAL-  AND  SHEAR-WAVE  PROFILES 
USED  IN  THE  APPLICATION  OF  THE 
SEMI AUTOMATED  SATELLITE  ANALYSIS  PROCEDURES 
TO  THE  32  DETECTED  INCLUSIONS 
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Figure  D-l.  Longitudinal-  and  Shear-Wave  Profiles  of  the  Inclusions 
in  Titanium  Specimen  No.  FML  97416  (Flaws  11  through  13) 
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Figure  D- 


.  Longitudinal-  and  Shear-Wave  Profiles  of  the  Inclusions 
in  Titanium  Specimen  No.  FML  97416  (Flaws  14  and  15) 
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Figure  D-3.  Longitudinal-Wave  Profiles  of  the  Inclusions  in 
IN100  Specimen  Nos.  FML  97474-1  through  -8 
(Flaws  16  through  23) 
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[■  Figure  D-5.  Longitudinal-Wave  Profiles  of  the  Inclusions 

K  •  in  IN100  Specimen  Nos.  FML  97474-1  through  -8 

>.-•  (Flaws  32  through  37) 
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Figure  D-6.  Longitudinal-Wave  Profiles  of  the  Inclusions 
in  IN100  Specimen  Nos.  FML  97474-1  through  -8 
(Flaws  38  through  40) 
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